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Glossary 

Bohr hypotheses: Set of hypotheses that Bohr introduced to explain the sta- 
bility OF THE atom. 

Coulomb potential: Potential that holds the electrons attached to the nu- 
cleus. 

Creation and Annihilation operators: Mathematical objects that create or an- 
nihilate PARTICLES. 

Dynamical phase: It is the factor that contains the time dependence of the 

WAVE function. 

Hydrogen atom: The simplest atom in nature, consisting of one electron re- 
volving around a proton. 

Laplace operator: A partial differential operator that contains the second 
partial derivatives with respect to the space coordinates. 

Magnetic quantum number: A quantum number that is associated to the direc- 
tion of projection of the of the angular momentum 

Orbital quantum number: The number associated to the quantization of the 
orbital motion 

Principal quantum number: The number associated to the radial solution of 
the Schrodinger equation and defines the energy of the allowed orbits. 

Quantized elec^tron orbits: The ele(;;trons can move in the atom only in very 

SPECIFIC orbits. 

Quantum Harmonic oscill.ator: The quantum analogous model for a potential 
WHOSE restoring force is proportional to the displacement. 

Quantum Mecilanic's: Theory of the laws that rule the mk;roscopic world, i.e. 
at the atomic sizes and less. 

Rutherford planetary atomic model: model to describe the atom in which the 

NUCLEUS is at THE CENTER AND THE ELECTRONS MOVE AROUND IT. 

Schrodinger equation: The differential equation that describes the movement 
OF atomic particles. 

Square well potential: A one-dimensional model of a potential with vertical 

WALLS 

Standing wave: The solutions to the Schrodinger equation that are stationary. 
Zeeman effect: The effect of splitting the electronic energy levels when an 
atom is immersed in a magnetic field. 



Summary 

In the early d.ays of the 20th century a set of important observations in atomic 
and molec;ular physic's could not be explained on the basis of the l.aws of classical 
physics. One of the main findings was the emission of light by exc;ited atoms with 

VERY particular FREQUENCIES. TO EXPLAIN THOSE FINDINGS A NEW DEVELOPMENT IN PHYSICS 
WAS NECJESSARY, NOW KNOWN AS QUANTUM MECHANICS. In PARTICULAR, THE CONCEPT OF 
STATIONARY STATES WAS INTRODUCED BY NiELS BOHR, IN 1913, IN ORDER TO EXPLAIN THOSE 
OBSERVATIONS AND THE STABILITY OF ATOMS. ACCORDING TO E.G. KeMBLE (1929), THE 
EXISTENCE OF DISCRETE ATOMIC AND MOLECULAR ENERGY LEVELS BROUGHT INTO MECHANICS 
A NEW KIND OF ATOMICITY SUPERPOSED ON THE ATOMICITY OF ELECTRONS AND PROTONS. We 
REVIEW HERE IN A HISTORICAL CONTEXT THE TOPIC OF STATIONARY STATES IN THE QUANTUM 
WORLD, INCLUDING THE GENERALIZATION TO THE PRIMARY IDEAS. We ALSO DISCUSS THE 

stationary states in one dimensional parabolic wells and the three dimensional 
Coulomb and parabolic cases. 
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1 Introduction 



At the beginning of the 20th century, some experimental observations in atomic and 

MOLECULAR PHYSICS WERE IMPOSSIBLE TO EXPLAIN ON THE BASES OF CLASSICAL PHYSICS. It 
WAS NECESSARY TO INTRODUCE REVOLUTIONARY CONCEPTS THAT LEAD TO THE FOUNDATION 
OF QUANTUM MECHANICS. In THIS CONTEXT THE CONCEPT OF STATIONARY STATES PLAYED AN 
ESSENTIAL ROLE IN THE DEVELOPMENT OF NEW IDEAS THAT STARTED TO EXPLAIN THE ATOMIC 
WORLD. 

In 1908 J.R. Rydberg and W. Ritz studied in detail the spectra of the light emit- 
ted BY EXCITED ATOMS. THEY FOUND THAT THE SPECTRA CONSISTED OF A SET OF DEFINED 

lines of particular wavelengths. furthermore, the set of spectroscopic lines were 
dependent only on the atom under study. through the so-called combination prin- 
ciple they put the data in a most systematic form. their principle states that the 
frequency of a particular spectral line can be expressed as a difference between 
some members of the set of frequency lines. 

These findings could not be explained by the accepted atomic model at that time, 
PROPOSED BY J.J. Thomson, claiming that the electrons were embedded in a positive 
charged cloud, whose extent was determined by the atomic radius. That model 
could not explain also the data obtained by H.W. Geiger and E. Mardsen, who 
under the supervision of Rutherford, were studying the interaction of charged 
o-particles with gold foils [t] . they observed that a considerable fraction of the 
a-particles was deflected by large angles. this effect could not be attributed 
to the electrons since they are much less massive. thus, they concluded that the 
source of deflection must be the positive charge concentrated in a much smaller 

VOLUME THAN THE ONE GENERATED BY THE ATOMIC RADIUS. In 1911 RUTHERFORD PROPOSED 

a new atomic model which assumed that all the positive charge is located at the 
center of the atom with a very dense distribution with a radius much smaller that 
the atomic one. the electrons then would circulate around the nucleus in a similar 
way as the planets move around the sun. 

Although Rutherford's planetary atomic model explained qualitatively well the 
deflection of a-particles, it had two major deficiencies. first it could not account 
for the spectra of radiation from atoms, which was not continuous but discrete. 
The other major problem was that, according to electrodynamics, an electron 
moving around the atom is under a constant acceleration, must radiate energy. 
This fact would lead to a situation in which the electron would loose energy 
continuously and would collapse with the nucleus. 

2 Stationary Orbits in Old Quantum Mechanics 
2.1 Quantized Planetary Atomic Model 

In 1911, THE TWO-AND-A HALF-THOUSAND-YEAR-OLD PHILOSOPHICAL CONCEPT OF ATOM TURNED 
INTO A SCIENTIFIC MATTER WHEN RUTHERFORD 'S PLANETARY ATOMIC MODEL EMERGED FROM 
THE INTERPRETATION OF THE EXPERIMENTAL DATA ON THE SCATTERING OF a PARTICLES [T|. 
The curious fact that has been noticed while these particles were SHUT TO GOLD FOILS 
WAS that some of them bounced AS IF THEY WERE COLLIDING WITH VERY MASSIVE OBJECTS. 
To EXPLAIN THESE FINDINGS RUTHERFORD PROPOSED THAT THE ATOM WAS COMPOSED BY A 
POSITIVE CENTRAL MASSIVE NUCLEUS AND THE ELECTRONS WERE REVOLVING AROUND IT, I.E. 
VERY SIMILAR TO A MINIATURE SOLAR SYSTEM. HOWEVER, THIS FAMOUS MODEL WAS NOT ELEC- 
TRODYNAMICALLY VIABLE. ATOMIC STABILITY WAS SIMPLY NOT ASSURED FOR RuTHERFORD'S 
SEMIEMPIRIC MODEL, SINCE ACCELERATED CHARGES RADIATE ENERGY AND THE ELECTRONS 
MOVING AROUND THE NUCLEUS WOULD LOSS ENERGY AND EVENTUALLY COLLAPSE WITH THE 
NUCLEUS. 

Another important set of empirical set of data, is that obtained from the emis- 
sion OF LIGHT BY EXCITED ATOMS. It WAS OBSERVED THAT THE LIGHT EMITTED HAD A VERY 
CHARACTERISTIC FREQUENCIES AND WAS A FOOTPRINT FOR EACH ATOM. THESE OBSERVATIONS 
WERE PUT IN A SYSTEMATIC FORM IN 1908 THROUGH THE SO-CALLED COMBINATION PRINCIPLE 
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FORMULATED BY J.R. RyDBERG AND W. RiTZ. THEIR PRINCIPLE SAYS THAT THE FREQUENCY 
OF A SPECTRAL EMISSION OR ABSORPTION LINE CAN BE EXPRESSED AS A DIFFERENCE BETWEEN 
THE MEMBERS OF A SET OF WELL DEFINED FREQUENCY TERMS. RuTHERFORD'S MODEL WAS 
COMPLETELY SILENT ON THE DYNAMICAL ORIGIN OF THE SPECTRAL LINES. It WAS THE GREAT 
MERIT OF Bohr to formulate in 1913 the HYPOTHESES, OR POSTULATES, THAT COULD AL- 
LOW THE EXPLANATION OF THE ATOMIC SPECTRAL LINES BASED ON THE PLANETARY ATOMIC 
STRUCTURE. 



2.2 Bohr's Hypotheses and Quantized Circular Orbits 

The HYPOTHESES THAT BOHR ADDED TO THE RUTHERFORD MODEL IN ORDER TO EXPLAIN THE 
SPECTROSCOPIC INFORMATION ARE THE FOLLOWING \2\ 

1. An ATOM CAN EXIST ONLY IN SPECIAL STATES WITH DISCRETE VALUES OF 
ENERGY. In other WORDS, THE ELECTRONS MOVING AROUND AN ATOM CAN 
BE FOUND ONLY IN CERTAIN SPECIAL ORBITS THAT BOHR CALLED stationary 
states. 

2. When an atom makes a transition from one stationary state to an- 
other, IT emits or absorbs radiation whose frequency u is given by 

THE frequency CONDITION 

hu = Ei- E2 , (1) 

where El AND E2 ARE THE ENERGIES OF TWO STATIONARY STATES. 

3. In the STATIONARY STATES, THE ELECTRONS MOVE ACCORDING TO THE 
LAWS OF CLASSICAL THEORY. HOWEVER, ONLY THOSE MOTIONS ARE PER- 
FORMED FOR WHICH THE FOLLOWING quantum condition IS FULFILLED 

pdq = nh, (n = 1, 2, 3, ) , (2) 

WHERE p IS THE MOMENTUM OF THE ELECTRON AND q IS ITS COORDI- 
NATE ALONG THE STATIONARY ORBIT. THE INTEGRATION SHOULD BE TAKEN 
ALONG THE ORBIT OVER ONE PERIOD OF THE CYCLIC MOTION. 

Bohr's theory claimed that those frequency terms, when multiplied by h, give 

DISTINCT ENERGY LEVELS IN WHICH THE ELECTRONS MOVE AROUND THE NUCLEUS. THIS MEANT 
THAT THESE WERE THE ONLY POSSIBLE STATES IN WHICH THE ELECTRONS IN THE ATOM COULD 
EXIST. 

Let us ASSUME that an electron in A Hydrogen atom is revolving around the 

NUCLEUS ON A CIRCULAR ORBIT ACCORDING TO THE NEWTONIAN EQUATIONS OF MOTION. FOR 
A CIRCULAR ORBIT, THE ABSOLUTE VALUE OF THE MOMENTUM p IS CONSTANT AND THEN THE 
QUANTUM HYPOTHESIS (3) LEADS TO 

p ■ 2-110 = nh, (n = 1,2,3, ...) (3) 

where a is the radius of the orbit. thus, a is given by the value of the momentum 
that can be obtained from the balance between the centrifugal force and the 
Coulomb force, i.e.. 



ma A-KtQa^ 
Combining the two equations, one obtains 



(4) 



enh'^n^ , , , 

a^ = - ^ (n= 1,2,3,...) . (5) 

The latter formula gives the radii of the quantized electron circles in the 
hydrogen atom. In particular, ai = an = ^° ■> , is known as the Bohr radius and is 

TAKEN AS AN ATOMIC LENGTH UNIT. 
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2.3 From Quantized Circles to Elliptical Orbits 



Wilson [3] and Sommerfeld [4] extended Bohr's ideas to a large variety of atomic 

SYSTEMS between 1915 AND 1916. 

The main idea is that the only classical orbits that are allowed AS STATIONARY 
STATES ARE THOSE FOR WHICH THE CONDITION 

<^Pkdqk=nkh k = l,...,n, (6) 

WITH rtk A POSITIVE INTEGER, IS FULFILLED. THE WEAK THEORETICAL POINT IS THAT 
IN GENERAL THESE INTEGRALS CAN BE CALCULATED ONLY FOR conditionally periodic systems, 
BECAUSE ONLY IN SUCH CASES A SET OF COORDINATES CAN BE FOUND, EACH OF WHICH GOES 
THROUGH A CYCLE AS A FUNCTION OF THE TIME, INDEPENDENTLY OF THE OTHERS. SOMETIMES 
THE COORDINATES CAN BE CHOSEN IN DIFFERENT WAYS, IN WHICH CASE THE SHAPES OF THE 
QUANTIZED ORBITS DEPEND ON THE CHOICE OF THE COORDINATE SYSTEM, BUT THE ENERGY 
VALUES DO NOT. 

In PARTICULAR, WHEN THE 3D POLAR COORDINATES ARE EMPLOYED, EQ. (6) GIVES THE 

Sommerfeld ellipses characterized by 

1^ Prdr = Urh , <j) pQdO = Ugh , <j> p^d4i = n^h . (7) 

Now, since p0 is a constant, one gets immediately the 'quantization' of the 
angular momentum of the ellipse along the z axis 

mh , , 

P0 = — , m = ±l,±2,---. (8) 

The quantum number m was called the magnetic quantum number by Sommer- 
feld WHO USED IT AS A MEASURE OF THE DIRECTION OF THE ORBIT WITH RESPECT TO THE 
MAGNETIC FIELD AND THUS EXPLAINING THE ZeEMAN EFFECT, I.E., THE SPLITTING OF THE 
SPECTROSCOPIC LINES IN A MAGNETIC FIELD. UNLESS FOR THE VALUE m = WHICH IS CON- 
SIDERED AS UNPHYSICAL, THIS 'OLD' m IS PRACTICALLY EQUIVALENT WITH SCHRODINGER'S m, 

which mathematically is the azimuthal separation constant but has a similar inter- 
pretation. 

Interestingly, and this is sometimes a source of confusion, the 'old' azimuthal 

QUANTUM NUMBER IS DENOTED BY k AND IS THE SUM OF ng AND m. It GIVES THE SHAPE OF THE 
ELLIPTIC ORBIT ACCORDING TO THE RELATIONSHIP | = 2 ^ WHERE n = Ur+k, ESTABLISHED BY 

Sommerfeld. Actually, this k is equivalent to Schrodinger's orbital number I plus 

1, BUT AGAIN THEIR MATHEMATICAL ORIGIN IS QUITE DIFFERENT. 



2.4 Experimental Proof of the Existence of Atomic Sta- 
tionary States 

The existence of discrete atomic energy levels was evidenced for THE FIRST TIME 

BY J. Franck and G. Hertz in 1914 [5]. They observed that when an electron col- 
lides WITH AN ATOM (MERCURY IN THEIR CASE), A TRANSFER OF A PARTICULAR AMOUNT OF 
ENERGY OCCURRED. THIS ENERGY TRANSFER WAS RECORDED SPECTROSCOPICALLY AND CON- 
FIRMED Bohr's hypotheses that atoms can absorb energy only in quantum portions. 
Even today, the experiment is preferentially done either with mercury or neon 
TUBES. From the spectroscopic evidence, it is known that the excited mercury va- 
por EMITS ultraviolet RADIATION WHOSE WAVELENGTH IS 2536 A, CORRESPONDING TO A 
photon ENERGY hv EQUAL TO 4.89 eV. 

The FAMOUS FrANCK-HeRTZ curves represent the ELECTRON CURRENT VERSUS THE 
ACCELERATING POTENTIAL, SHOWN IN FiG. 3. THE CURRENT SHOWS A SERIES OF EQUALLY 
SPACED MAXIMA (AND MINIMA) AT THE DISTANCE OF ~ 4.9 V. THE FIRST DIP CORRESPONDS 
TO ELECTRONS THAT LOSE ALL THEIR KINETIC ENERGY AFTER ONE INELASTIC COLLISION WITH 
A MERCURY ATOM, WHICH IS THEN PROMOTED TO ITS FIRST EXCITED STATE. THE SECOND DIP 
CORRESPONDS TO THOSE ELECTRONS THAT HAVE THE DOUBLE AMOUNT OF KINETIC ENERGY AND 
LOSES IT THROUGH TWO INELASTIC COLLISIONS WITH TWO MERCURY ATOMS, AND SO ON. AlL 
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Figurc 1: Bohr-Sommerfeld electron orbits for n = 1, 2, and 3, and the allowed 
values for I. 
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Figure 2: Spatial quantization of Bohr-Sommerfeld orbits for azimuthal numbers 
k = 1,2, and 3. 
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THESE EXCITED ATOMS EMIT THE SAME RADIATION AT ~ 2536 A. BUT WHICH IS THE 'FIRST' 
EXCITED STATE OF MERCURY? It IS SPECTROSCOPIOALLY DENOTED BY Pi IN FiG. ((JJl. NOTICE 
THAT THE OTHER TWO P STATES CANNOT DECAY TO THE GROUND STATE ^ Sq BECAUSE THE 

dipole emission is forbidden for them and therefore they are termed metastable. 
More details, such that the observed peak separation depends on the geometry of 

THE TUBE AND THE HG VAPOR PRESSURE, ARE EXPLAINED IN THE READABLE PAPER OF HANNE 



3 Stationary States in Wave Mechanics 
3.1 The Schrodinger Equation 

According to L. Pauling and E. Bricjht Wilson Jr. [7], already in the years 1920- 
1925 A decline of the 'old quantum theory' as the Bohr-Sommerfeld atomic theory 
IS historically known and which is based on the 'whole number' quantization of 

cyclic orbits was PATENT: ONLY VERY RECENTLY THERE IS SOME REVIVAL, ESPECIALLY IN 
THE MOLECULAR CONTEXT [§]. BUT IN 1925, A QUANTUM MECHANICS BASED ON THE MATRIX 
CALCULUS WAS DEVELOPED BY W. HeISENBERG, M. BORN, AND P. JORDAN AND THE BEST WAS 
TO COME IN 1926 WHEN SCHRODINGER IN A SERIES OF FOUR PAPERS DEVELOPED THE MOST 
EMPLOYED FORM OF QUANTUM MECHANICS, KNOWN AS wave mechanics. THE ADVANTAGE OF 
HIS THEORY OF ATOMIC MOTION IS THAT IT IS BASED ON STANDARD (pARTIAl) DIFFERENTIAL 
EQUATIONS, MORE EXACTLY ON THE StURM-LiOUVILLE THEORY OF SELF-ADJOINT LINEAR DIF- 
FERENTIAL OPERATORS. Schrodinger starts the first paper in the 1926 series with the 

FOLLOWING sentence 0: 

In THIS PAPER I wish TO CONSIDER, FIRST THE SIMPLEST CASE OF THE HYDROGEN 

atom, and show that the customary quantum conditions can be replaced 
by another postulate, in which the notion of 'whole numbers', merely 
as such, is not introduced. 

Indeed, he could obtain the basic equation of motion in nonrelativistic quantum 
mechanics, the so called schrodinger equation for the wavefunctions <e'(a;, t), and 
provided several analytic applications, among which was the hydrogen atom. the 

ORIGINAL DERIVATION IS BASED ON THE VARIATIONAL CALCULUS WITHIN THE StURM-LiOUVILLE 
APPROACH AND WAS GIVEN EIGHTY YEARS AGO. In HIS FIRST PAPER OF 1926, SCHRODINGER 
STATES THAT THE WAVEFUNCTIONS SHOULD BE SUCH AS TO MAKE THE 'HAMILTON INTEGRAL' 

JsM = J {fi^T{q, d-4>/dq) + (It , (9) 

STATIONARY SUBJECT TO THE NORMALIZING CONDITION / ^^dr = 1 WHICH CAN BE INCORPO- 
RATED THROUGH THE LAGRANGE MULTIPLIERS METHOD. THE EuLER-LaGRANGE EQUATION OF 
THE FUNCTIONAL Js ['I'] IS THE TIME-DEPENDENT SCHRODINGER EQUATION 

ih— = H<f . (10) 
When the wave function of the time-dependent Schrodinger equation is written in 

THE multiplicative FORM '^(x, t) = 1p(x)J^(t) ONE OBTAINS A COMPLETE SEPARATION OF THE 
SPACE AND TIME BEHAVIORS OF ^: ON ONE SIDE, ONE GETS THE STATIONARY SCHRODINGER 
EQUATION FOR 1p{x), 

- —-^ + V{x)i, = E^ , (11) 

ZiTTh CLX 

AND ON THE OTHER SIDE, THE SIMPLE TIME-DEPENDENT EQUATION FOR THE LOGDERIVATIVE OF 

in''-}^ = E (12) 

dt 

This decoupling of space and time components is possible whenever the potential 

ENERGY IS independent OF TIME. 

The space component has the form of a standing-wave equation. Thus, it is 
correct to regard the time-independent schrodinger equation as a wave equation 
from the point of view of the spatial phenomenology. 
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Figure 3: (a) Schematic diagram of the Franck-Hertz experiment, where the tube is 
filled with a gas of Mercury ; (b) typical plot recorded in a Franck-Hertz experiment 
with mercury, showing the periodic maxima. 
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Figure 4: Simplified Mercury's level diagram in the low-energy region in which the 
so-called hyperfine structure is neglected. The numbers 6 and 7 are Bohr's 'whole 
numbers' or Schrodinger's principal quantum numbers. 

3.2 The Dynamical Phase 

Furthermore, the time-dependence is multiplicative and reduces to a modulation 
of the phase of the spatial wave given by 

T = e-'-E;*/^ = cos{Et/h) - isin{Et/h) . (13) 

The PHASE FACTOR = e^'-^'/'' IS KNOWN AS THE dynamical phase. In recent times, 

OTHER PARAMETRIC PHASES HAVE BEEN RECOGNIZED TO OCCUR, E.G., THE BeRRY PHASE. THE 
DYNAMK;AL PHASE IS A harmonic; OSCILLATION WITH ANCiULAR FRECJUENCY w = E/h AND 
PERIOD T = h/E. In OTHER WORDS, A SCMRODINCiER WAVEFUNCTION IS FLICKERING FROM 
POSITIVE THROUGH IMAGINARY TO NECJATIVE AMPLITUDES WITH A FRECJUENCY PRGPORTIONiVL 
TO THE ENERGY. AlTHOUCSH IT IS A WWE OF CONSTANT ENERCSY IT IS NOT STATIONARY BE- 
CAUSE ITS PHASE IS TIME DEPENDENT (PERIODIC). HCJWEVER, A REMARKABLE FAC;T IS THAT 
THE PRODUCT I.E., THE MODULUS OF THE SCHRODINGER CONSTANT-ENERGY WAVES 

REMAINS CONSTANT IN TIME 

**^f = (14) 

It IS IN THE SENSE OF THEIR CONSTANT MODULUS THAT SCHRODINGER CONSTANT-ENERGY 
WAVES ARE CALLED STATIONARY STATES. 

3.3 The Schrodinger Wave Stationarity 

Thus, non-relativistic quantum stationarity refers tc3 wwes of c;onstant enerc;y 
and constant modulus, but not of constant phase, which can occur as solutions of 
schrodinger equation for time-independent potentials. in the schrodinger frame- 
work, the dynamical systems are usually assumed to exist in stationary states (or 

WAVES OF THIS TYPE). It IS WORTH NOTING THAT THE PREFERRED TERMINOLOGY IS THAT 
OF STATES AND NOT OF WAVES. THIS IS DUE TO THE FACT THAT BEING OF CONSTANT EN- 
ERGY THE SCHRODINGER STATIONARY WAVES DESCRIBE PHYSICAL SYSTEMS IN CONFIGURATIONS 
(or states) OF CONSTANT ENERGY WHICH CAN THEREFORE BE NATURALLY ASSOCIATED TO THE 
TRADITIONAL CONSERVATIVE HAMILTONIAN SYSTEMS. MOREOVER, THE LOCALIZATION OF THESE 
WAVES CAN BE ACHIEVED BY IMPOSING APPROPRIATE BOUNDARY CONDITIONS. 
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3.4 Stationary Schrodinger States and Classical Orbits 



In the Schrodinger theory, a single stationary state does not correspond to a 

CLASSICAL orbit. ThIS IS WHERE THE SCHRODINGER ENERGY WAVES DIFFER THE MOST FROM 

Bohr's theory which is based on quantized classical cyclic trajectories. To build 
A wave entity closer to the concept of a classical orbit, one should use superposi- 
tions OF MANY stationary STATES, INCLUDING THEIR TIME DEPENDENCE, I.E., WHAT IS KNOWN 
AS WAVE PACKETS. OnLY MONOCHROMATIC PLANE WAVES OF ANGULAR FREQUENCY U CORRE- 
SPONDS THROUGH THE BASIC FORMULA E = hjj TO A WELL-DEFINED ENERGY E OF THE 'CLAS- 
SICAL' PARTICLE BUT UNFORTUNATELY THERE IS NO RELATIONSHIP BETWEEN THE WAVEVECTOR 
k AND THE MOMENTUM p OF THE CORRESPONDING PARTICLE SINCE A PLANE WAVE MEANS ONLY 
THE PROPAGATION AT CONSTANT (pHASe) VELOCITY OF INFINITE PLANES OF EQUAL PHASE. In 
OTHER WORDS, A CRITERIUM FOR LOCALIZATION IS REQUIRED IN ORDER TO DEFINE A CLASSICAL 
PARTICLE BY MEANS OF A WAVE APPROACH. 

In THE ONE-DIMENSIONAL CASE, A WAVE PACKET IS CONSTRUCTED AS FOLLOWS 

/ + 00 
/{fc')e'(*^ *)dfc' , (15) 

WITH OBVIOUS GENERALIZATION TO MORE DIMENSIONS. If f{k') IS WRITTEN IN THE POLAR 

FORM F{k')e'-°' AND F IS CHOSEN WITH A PRONOUNCED PEAK IN A WAVENUMBER REGION OF 

EXTENSION Ak AROUND THE POINT k' = k, THEN THE WAVE PACKET IS LOCALIZED IN A SPATIAL 

REGION OF EXTENSION A.X ft! ^jT SURROUNDING THE "CENTER OF THE WAVEPACKET" . THE 

LATTER IS EQUIVALENT TO THE CONCEPT OF MATERIAL POINT IN CLASSICAL MECHANICS AND 

TRAVELS UNIFORMLY WITH THE GROUP VELOCITY D„ = . THIS IS THE VELOCITY THAT CAN 

y dk 

BE IDENTIFIED WITH THE PARTICLE VELOCITY D = 4^ IN CLASSICAL MECHANICS AND WHICH 
LEADS TO THE DE BrOGLIE FORMULA p = Hk = j. 



3.5 Stationary States as Sturm-Liouville Eigenfunctions 

The mathematical basis of Schrodinger wave mechanics is the Sturm-Liouville (SL) 

THEORY of self-adjoint LINEAR DIFFERENTIAL OPERATORS ESTABLISHED IN THE 19TH CEN- 
TURY, MORE SPECIFICALLY THE SL EIGENVALUE PROBLEM, WHICH IS TO FIND SOLUTIONS OF A 
DIFFERENTIAL EQUATION OF THE FORM 

d 



CUn{x) _ 

dx 



dUn{x) 

p{x) 



dx 



■ q(x)Un(x) = -flnW{x)Un{x) (16) 



SUBJECT TO SPECIFIED BOUNDARY CONDITIONS AT THE BEGINNING AND END OF SOME INTERVAL 
(a, b) ON THE REAL LINE. E.G. 

Au{a) + Bu' (a) = (17) 
Cu{b) + Du' (b) = (18) 

WHERE A, B, C, D ARE GIVEN CONSTANTS. THE DIFFERENTIAL OPERATOR C IN (I16I I IS RATHER 
GENERAL SINCE ANY SECOND ORDER LINEAR DIFFERENTIAL OPERATOR CAN BE PUT IN THIS FORM 
AFTER MULTIPLICATION BY A SUITABLE FACTOR. THE BOUNDARY CONDITIONS ARE ALSO RATHER 
GENERAL INCLUDING THE WELL-KNOWN DiRICHLET AND NEUMANN BOUNDARY CONDITIONS AS 
PARTICULAR CASES BUT OTHER POSSIBILITIES SUCH AS PERIODIC BOUNDARY CONDITIONS 

u{x) = u{x + b- a) (19) 

COULD BE OF INTEREST IN SOME CASES, ESPECIALLY' FOR ANGULAR VARIABLES. 

The SL eigenvalue problem is an infinite dimensional generalization of the FINITE 
DIMENSIONAL MATRIX EIGENVALUE PROBLEM 

Mu = ^lu (20) 

WITH M AN n X n MATRIX AND 1i AN n DIMENSIONAL COLUMN VECTOR. AS IN THE MATRIX 
CASE, THE SL EIGENVALUE PROBLEM WILL HAVE SOLUTIONS ONLY FOR CERTAIN VALUES OF THE 
EIGENVALUE Pn- THE SOLUTIONS U„ CORRESPONDING TO THESE p,„ ARE THE EIGENVECTORS. 
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For the finite dimensional case with an n x n matrix M there can be at most n 

LINEARLY INDEPENDENT EIGENVECTORS. FOR THE SL CASE THERE WILL IN GENERAL BE AN 
infinite SET OF EIGENVALUES /i„ WITH CORRESPONDING EIGENFUNCTIONS Un{x) . 

The DIFFERENTIAL EQUATIONS DERIVED BY SEPARATING VARIABLES ARE IN GENERAL OF THE 
SL FORM, THE SEPARATION CONSTANTS BEING THE EIGENVALUE PARAMETERS fl. THE BOUNDARY 
CONDITIONS II17I18I I ARE DETERMINED BY THE PHYSICAL APPLICATION UNDER STUDY. 

The solutions Un{x), fin OF A SL EIGENVALUE PROBLEM HAVE SOME GENERAL PROPERTIES 
OF BASIC IMPORTANCE IN WAVE (QUANTUM) MECHANICS. 

If u{x) and v(x) are arbitrary twice differentiable solutions of a SL operator, 

THEN by integrating BY PARTS 

'' du dv^^ , . 

(21) 



dx[vCu — uCv] = p j V u — 

dx dx 



An operator C which satisfies II21II is said to be self-adjoint. Any second order 

LINEAR DIFFERENTIAL OPERATOR CAN BE PUT IN THIS SELF ADJOINT FORM BY MULTIPLICATION 
BY A SUITABLE FACTOR. 

It IS EASY TO SHOW THAT FOR FUNCTIONS u{x), v{x) SATISFYING BOUNDARY CONDITIONS 
OF THE STANDARD SL FORM OR THE PERIODIC BOUNDARY CONDITION II19II THE RIGHT HAND 
SIDE OF II21II VANISHES. FOR BOTH OF THIS CASE WE THEN HAVE 

f-b rh 

I dxvCu = / dxuCv . (22) 

•J a J a 

Consider now two different eigenfunctions Un{x), Um{x) belonging to different 

EIGENVALUES A„ 7^ : 

CUn{x) = —fl„w{x)Un{x) , (23) 
CUm(x) = - fJLmW{x)Um{x) . (24) 

Multiplying the first equation by Um and the second one by «„, integrating and 

SUBTRACTING, WE FIND: 

/h rh 
dx[Um.CUn - UnLUm] = -{fJ-n - flm) / dxw {x)Un {x)Um (x) . (25) 
•J a 

The LEFT hand side will vanish for EITHER SET OF BOUNDARY CONDITIONS WE CONSIDER 
HERE, SO FOR EITHER OF THESE CASES. WE FIND THE ORTHOGONALITY CONDITION 



b 

dxw{x)Un{x)Um{x) = , fJ.n ^ flm ■ (26) 



Two FUNCTIONS U„{x), Um{x) SATISFYING THIS CONDITION ARE SAID TO BE ORTHOGONAL 
WITH WEIGHT w{x). MOREOVER, IF w{x) IS NON NEGATIVE, ONE CAN INTRODUCE THE SL 
NORMALIZATION OF THE Un AS FOLLOWS: 



b 

dxw{x)[unix)f = 1 . (27) 



The most important property of the eigenfunctions of a SL problem is that they 

FORM a complete SET. THIS MEANS THAT AN ARBITRARY FUNCTION tp{x) CAN BE EXPANDED 
IN AN INFINITE SERIES OF THE FORM 

i>{x) = ^a,i?i„(x) . (28) 
n 

The expansion coefficients a„ are determined by multiplying II28I I by w{x)un{x), in- 
tegrating TERM BY TERM, AND USING THE ORTHOGONALITY RELATION 

_ dxw{x)Un(x)i}(x) ^^^^ 



Ja d.Xw(x)[Un{x)] 



2 
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According to Courant and Hilbert [To], every piecewise continuous function 
defined in some domain with a square-integrable first derivative may be expanded in 
an eigenfunction series which converges absolutely and uniformly in all subdomains 
free of points of discontinuity; at the points of discontinuity it represents (like the 
Fourier series) the arithmetic mean of the right and left hand limits. (This theorem 

DOES not require THAT THE FUNCTIONS EXPANDED SATISFY THE BOUNDARY CONDITIONS.) 



4 The Infinite Square Well: The Stationary States 
Most Resembling the Standing Waves on a 
String 

We HAVE already COMMENTED THAT THE TIME-INDEPENDENT SCHRODINGER EQUATION HAS 
THE FORM OF A STANDING- WAVE EQUATION. THIS IS A VERY INSTRUCTIVE ANALOGY AND AL- 
LOWS TO OBTAIN THE CORRECT ENERGY VALUES FOR THE CASE OF THE INFINITE SQUARE WELL 
USING ONLY THE DE BrOGLIE WAVE CONCEPT WITHOUT EVEN INTRODUCING THE SCHRODINGER 
EQUATION nil . 

We REMIND THE TREATMENT OF THE STRING STANDING WAVES IN THE CASE OF A FINITE 
HOMOGENEOUS STRING OF TOTAL LENGTH L. PROVIDED THAT THE ORIGIN OF THE COORDINATE 
SYSTEM IS PLACED IN ITS CENTER AND THE X DIRECTION IS CHOSEN PAR ALLEL TO IT, THE SPACE 
DEPENDENT PART OF ITS STANDING WAVES IS GIVEN BY 

u{x) = A COS kx + B sin kx (30) 

WHERE X e ( — L/2, L/2) AND A, B, AND k ARE CONSTANTS. IMPOSING THE USUAL (DiRICHLET) 
BOUNDARY CONDITIONS 

n(-L/2) = u{L/2) = , (31) 

EQUATION l|30|l TAKES THE FORM 

, , f Acosfc„a; , IF n > AND ODD 

Unix) = < n ■ 1 r, (32) 

^ ' I BsmknX , IF n > AND EVEN ' 



AND 



X G (-L/2, L/2) (33) 

fcn = ^ , (n = 1,2,3,...) . (34) 
These functions - the normal modes of the string under consideration - form a 

COMPLETE set WITH RESPECT TO PHYSICALLY REASONABLE FUNCTIONS DEFINED WITHIN THE 
INTERVAL (I33II AND SATISFYING EQUATION (I31I I. I.E., ANY SUCH FUNCTION U (x) CAN BE WRITTEN 

AS THE Fourier series 

oo 

= CnU„{x) . (35) 

n=l 

We are READY NOW TO STATE THE ANALOGY, WHICH IS BASED ON THE FOLLOWING STEPS. 

• A STANDING DE BrOGLIE WAVE CORRESPONDS TO A QUANTUM PARTICLE STRICTLY CON- 
FINED TO THE REGION —L/2 < X < L/2, I.E., IN AN INFINITE SQUARE- WELL POTENTIAL 

SHOWN IN (Fig. ^ 

n-) = (°' "ll^f-^' (36) 

^ ' [ -oo , IF \x\ > L/2 . ^ ' 

Because of the form of this potential, it is assumed that there is no asymp- 
totic TAIL OF THE WAVE FUNCTIONS IN THE OUTSIDE REGIONS TO THE WELL, I.E., 
1p{x) = FOR \x\ > L/2. In physical terms, that means that the particle is COM- 
PLETELY LOCALIZED WITHIN THE WALL AND THE BOUNDARY CONDITIONS ARE OF THE 
DiRICHLET TYPE. 
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V(x) 



-L/2 



L/2 



Figure 5: The infinite square well potential. This case is mathematically the 
most analogous to the classical string standing waves. 



In the interior region, the SCHRODINGER equation is the simplest POSSIBLE: 

WHOSE SOLUTIONS (ENERGY EIGENFUNCTIONS) ARE: 

(2/L)l/2 sin(n7rx/L) , nEVEN , \x\ < L/2 , 
<t>„{x) = { (2/L)i/2cos(n7rx/L) , n ODD , \x\<L/2, (38) 



, |x| > L/2 , 

AND THE ENERGY EIGENVALUES ARE: 

2m 2mL2 

• The AMPLITUDE OF THIS STANDING WAVE IN THE nTH STATIONARY STATE IS PROPOR- 
TIONAL TO sm(n7rz/L). This corresponds strictly to the analogy with standing 

WAVES on a classical STRING. 

• The nTH standing wave is presented as a SUPERPOSITION OF TWO RUNNING WAVES. 

The wave travelling right, with wavelength A„ = 2L/n has de Broglie momen- 
tum fmn/L AND THE LEFT TRAVELLING WAVE HAS OPPOSED DE BrOGLIE MOMENTUM 

— fmn/L. The resulting energy is then quantized and given by 

Because of the extremely strong confinement of the infinite square well it 
seems that this case is only of academic interest. However, two-dimensional strong 
confinement of electrons by rings of adatoms (corrals) have been reported in the 
literature [12]. 
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Figure 6: The harmonic (parabohc) oscillator potential. 

5 ID Parabolic Well: The Stationary States of 
the Quantum Harmonic Oscillator 

5.1 The Solution of the Schrodinger Equation 

The harmonic oscillator (HO) is one of the fundamental paradigms of Physics. Its 

UTILITY resides IN ITS SIMPLICITY WHICH IS MANIFEST IN MANY AREAS FROM CLASSICAL PHYSIOS 
TO QUANTUM ELECTRODYNAMICS AND THEORIES OF GRAVITATIONAL COLLAPSE. 
It IS WELL KNOWN THAT WITHIN CLASSICAL MECHANICS MANY COMPLICATED POTENTIALS CAN 
BE WELL APPROXIMATED CLOSE TO THEIR EQUILIBRIUM POSITIONS a; BY HO POTENTIALS AS 
FOLLOWS 

(41) 



1 



Vix) ~ -V" {ai)ix - a,y 

For THIS CASE, THE CLASSICAL HAMILTONIAN FUNCTION OF A PARTICLE OF MASS m, OS- 
CILLATING AT THE FREQUENCY UJ HAS THE FOLLOWING FORM: 

^2 



H : 



1 



2_2 



1 — mu) X 

2 m 2 



(42) 



AND THE QUANTUM HAMILTONIAN CORRESPONDING TO THE CONFIGURATIONAL SPACE IS GIVEN 
BY 



H : 



(43) 



1/ dYl 22 d'^ 1 22 

—in — H — muj x = H — mo; x . 

2m \ dx J 2 2m dx^ 2 

Since we consider a time-independent potential, the eigenfunctions i/in and the 

EIGENVALUES ARE OBTAINED BY MEANS OF THE TIME-INDEPENDENT SCHRODINGER EQUATION 



Hljj = Elp 



(44) 
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For the ho Hamiltonian, the Schrodinger, equation is 

2mE m?u 



Defining the parameters 



2mE 



= 



A : 



muj 



THE Schrodinger equation becomes 

WHICH IS KNOWN AS WeBER'S DIFFERENTIAL EQUATION IN MATHEMATICS. TO SOLVE THIS 
EQUATION ONE MAKES USE OF THE FOLLOWING VARIABLE TRANSFORMATION 



A^x^lif = 



(45) 
(46) 
(47) 



y = \x 



(48) 



By CHANGING THE INDEPENDENT VARIABLE FROM X TO y, THE DIFFERENTIAL OPERATORS 
AND D? TAKE THE FORM 



d 



Dt 



(49) 



dx \dx J " ' ^ dx^ dx \dx dy J \dx'^ J \dx J ^ 

By applying these rules to the proposed TRANSFORMATION WE OBTAIN THE FOLLOW- 
ING DIFFERENTIAL EQUATION IN THE y VARIABLE 



y 



dy'^ 



1 dtp 



fc2 



1 



4A 4^ 



tp = 



AND, BY DEFINING: 



fc2 _ _B 
2\~ hu! 



WE GET AFTER DIVIDING BY y (i.E., y ^ 0) 



1 djjj 
2y dy 



2V 



i/) = 



(50) 
(51) 

(52) 



Let us TRY TO SOLVE THIS EQUATION BY FIRST DOING ITS ASYMPTOTIC ANALYSIS IN THE LIMIT 
y ^ OO IN WHICH THE EQUATION BEHAVES AS FOLLOWS 



(53) 



dy2 4"^ 
This equation has as solution 

V'oo(?;) = Aexp(|) +Bcxp(-|) . (54) 

The first term diverges in the limit y ^ oo. Thus we take A = 0, and keep only 
the attenuated exponential. We can now suggest that ■0 has the following form 



Plugging it in the differential equation for y ( Eq. (I52I i1 one gets: 
d'^ip (\ \ dip (k \\ 

The latter equation is of the confluent (Kummer) hypergeometric form 



WHOSE GENERAL SOLUTION IS 



!{z) = A iFi(a, c;2) + B ^l-=iFi(a- 



■l,2-c;2) 



(55) 
(56) 

(57) 
(58) 
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WHERE THE CONFLUENT HYPERGEOMETRIC FUNCTION IS DEFINED AS FOLLOWS 



iFi{a,c;z) ■ 



(59) 



By DIRECT COMPARISON OF EQS. II57I I AND II56II . ONE CAN SEE THAT THE GENERAL SOLUTION 
OF THE LATTER ONE IS 



■ A iFi{ a, -■,y 



B j/2 iFi 



1 3 
2' 2' 



(60) 
(61) 



K V 

If we keep THESE SOLUTIONS IN THEIR PRESENT FORM, THE NORMALIZATION CONDITION IS NOT 
SATISFIED BY THE WAVEFUNCTION. INDEED, BECAUSE THE \y\ OO ASYMPTOTIC BEHAVIOUR OF 

THE CONFLUENT HYPERGEOMETRIC FUNCTION IS iFi{a,C;y) rrl-n.l ^~^'^'^y~°' + TJ^^^y°'~'^ ' 



r(c-a) 

IT FOLLOWS FROM THE DOMINANT EXPONENTIAL BEHAVIOR THAT: 



const. e^2'y^^ 



(62) 



This leads to a dfvergence in the normalization integral, which physically is not 
acceptable. what one does in this case is to impose the termination condition for 
the series. that is, the series is cut to a finite number of n terms and therefore 
turns into a polynomial of order n. the truncation condition of the confluent 

HYPERGEOMETRIC SERIES iFi{a,C;z) IS a = —n, WHERE n IS A NONNEGATIVE INTEGER (I.E., 
ZERO IS included). 

We thus notice that asking for a finite normalization CONSTANT, (as already KNOWN, 
A necessary condition for THE PHYSICAL INTERPRETATION IN TERMS OF PROBABILITIES), 
leads us TO THE TRUNCATION OF THE SERIES, WHICH SIMULTANEOUSLY GENERATES THE QUAN- 
TIZATION OF ENERGY. 

In THE FOLLOWING WE CONSIDER THE TWO POSSIBLE CASES: 
1) a = —n AND B = 



The eigenfunctions are given by 

■0„(a;) = Nn cxp 

and THE energy IS: 



1 

~ 4 



hjj\2n + ■ 



■,\x' 



2) 



and a = 



= n - 



The EIGENFUNCTIONS ARE NOW 

V'n(a;) = Nn exp 



WHEREAS THE STATIONARY ENERGIES ARE 



En = fujJ 



(2n- 



: \x' 



(63) 
(64) 
(65) 

(66) 
(67) 
(68) 



The polynomials obtained by this truncation of the confluent hypergeometric se- 
ries ARE CALLED HeRMITE POLYNOMIALS AND IN HYPERGEOMETRIC NOTATION THEY ARE 

1 



H2n{n) = (-1)"^ iFl 



H.„_.(.) = , ,F. 



(69) 
(70) 
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We can now combine the obtained results (because some of them give us the 

EVEN cases and THE OTHERS THE ODD ONES) IN A SINGLE EXPRESSION FOR THE OSCILLATOR 
EIGENVALUES AND EIGENFUNCTIONS 

4)n{x) = A/'nexpf —\H„{y/Xx) (71) 

E„ = (^n + ^^ hw n = 0, 1, 2 . . . (72) 

The ho energy spectrum is equidistant, i.e., there is the same energy difference 
fitj between any consecutive neighbor levels. Another remark refers to the min- 
imum VALUE OF THE ENERGY OF THE OSCILLATOR; SOMEWHAT SURPRISINGLY IT IS NOT ZERO. 

This is considered by many people to be a pure quantum result because it is zero 

WHEN ?i — > 0. Eo = ^huj IS known AS THE ZERO POINT ENERGY AND THE FACT THAT IT IS 
nonzero is THE MAIN CHARACTERISTIC OF ALL CONFINING POTENTIALS. 



5.2 The Normalization Constant 

The normalization constant is usually calculated in the following way. The Her- 

MITE GENERATING FUNCTION e^t~* +2tx) jg MULTIPLIED BY ITSELF AND THEN BY e~^^ : 



->.a^\x(-s'+2s.)^x(-t'+2t.) ^ ^ e-^-^f^„(yAz)/^„(VAx•) 

m.n—O 



(73) 



Integrating over x' = \/~Xx on the whole real line, the cross terms of the double 
SUM drop out because of the orthogonality property 



2^ („\\2 I 



e~^-[H„{x')fdx' 



f 

J — < 



2"A"(st)'' 



WHERE THE PROPERTIES OF THE EULER GAMMA FUNCTION HAVE BEEN USED 

/■oo 

r{z) = 2 t^^-^di , Re(z) > 

Jo 

AS WELL AS THE PARTICULAR CASE Z = 1/2 WHEN r(l/2) = y'TT. 

By equating coefficients of like POWERS OF st IN (I74p . WE OBTAIN 



/oo 2 
e-^ [H„{x')]^dx' =2"n^^^n\ 
- oo 



This leads to 



X 1 

TT 2"n! 



(74) 



(75) 



(76) 



(77) 



5.3 Final Formulas for the HO Stationary States 

Thus, one gets the following normalized eigenfunctions (stationary states) of the 
one-dimensional harmonic oscillator operator 



X 1 

TT 2"n! 



cxp 



-Ax" 



H^iVXx) 



(78) 
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If the dynamical phase factor T is included, the harmonic oscillator eigenfunctions 

TAKES the following FINAL FORM 



A 1 
TT 2"n! 



exp 



[~i (ri + i j c^t - ^ j //„(v^x) . (79) 



5.4 The Algebraic Approach: Creation and Annihilation 
Operators and a 

There is another approach to deal with the HO besides the conventional one of 
solving the Schrodinger equation. It is the algebraic method, also known as the 

METHOD of creation AND ANNIHILATION (lADDER) OPERATORS. THIS IS A VERY EFFICIENT 
procedure, which CAN BE SUCCESSFULLY APPLIED TO MANY QUANTUM-MECHANICAL PROB- 
LEMS, ESPECIALLY WHEN DEALING WITH DISCRETE SPECTRA. 
Let us DEFINE TWO NONHERMITIC OPERATORS a AND : 



mu) I ip 
2h \ mui 



(80) 
(81) 



J. / muj / ip 

a' = , I X 

\ 2h \ moj 

These operators are known as annihilation operator and creation operator, 
respectively (the reason of this terminology will be seen in the following. 

Let US CALCULATE THE COMMUTATOR OF THESE OPERATORS 

[a,at] = — x + -i-,x--i- =—{-i[x,p]+i[p,x]) = l, (82) 

WHERE WE HAVE USED THE COMMUTATOR [x,p] = iH. THEREFORE THE ANNIHILATION AND 
CREATION OPERATORS DO NOT COMMUTE, SINCE WE HAVE [a, at] = 1. LET US ALSO INTRODUCE 
THE VERY IMPORTANT NUMBER OPERATOR N: 

N = a)a . (83) 

This operator is hermitic as one can readily prove using (ABy = B'^ A'' : 

Aft = (ata)"!" = at(a'l')t = ata = TV . (84) 

Considering now that 



H 1 

2h V ' m2w2 J 2ft" " ' "hi) 2 



^' + 1^1 +-[^.P] = — ' (85) 



we notice that the Hamiltonian can be written in a quite simple form as a linear 
function of the number operator 

H = nuj(^N . (86) 

The number operator bears this name because its eigenvalues are precisely the 
subindexes of the eigenfunctions on which it acts 

7V|n> = n\n) , (87) 

WHERE WE HAVE USED DiRAC'S KET NOTATION Vn = I"). 

Applying the number-form of the HO Hamiltonian in (I86I I to this ket, one gets 

H\n) = h^(^+]^ \n) , (88) 
WHICH DIRECTLY SHOWS THAT THE ENERGY EIGENVALUES ARE GIVEN BY 



h^[n+l] . (89) 
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Thus, this basic result is obtained through purely algebraic means. 

It is possible to consider the ket at|n} as an eigenket of that number opera- 
tor FOR WHICH THE EIGENVALUE IS RAISED BY ONE UNIT. In PHYSICAL TERMS, THIS MEANS 
THAT AN ENERGY QUANTA HAS BEEN PRODUCED BY THE ACTION OF at ON THE KET \n) . TlIIS 
ALREADY EXPLAINS THE NAME OF CREATION OPERATOR. SIMILAR COMMENTS WITH CJORRE- 
SPONDINC; CONCLUSION CAN BE INFERRED FOR THE OPERATOR a EXPLAINING THE NAME OF 
ANNIHILATION OPERATOR (AN ENERGY QUANTA IS ELIMINATED FROM THE SYSTEM WHEN THIS 

operator is put into action). 
Consequently, we have 

at|n) = Vn+T|n + 1) . (90) 

For the annihilation operator, following the same procedure one can get the 
following relation 

a\n) = Vn|n - 1> . (91) 
Let us show now that the values of n should be nonnegative integers. For this, 

we EMPLOY THE POSITIVITY REQUIREMENT FOR THE NORM OF THE STATE VECTOR a\n) . ThE 
latter CONDITION TELLS US THAT THE INNER PRODUCT OF THE VECTOR WITH ITS ADJOINT 

(a|n))t (= (n|at) should always be nonnegative 

((n|at) . (a|n)) > . (92) 

This relationship is in fact the expectation value of the number operator 

(n|ata|n) = (n|iV"|n) = (n|a+-v/n|n - 1) = {n\VnVn\n) = n>0 . (93) 

Thus, n cannot be negative. It should be a positive integer since, if that would not 
be the case, by applying iteratively the annihilation operator a sufficient number 
of times we would be led to imaginary and negative eigenvalues, which would be a 
contradiction to the positivity of the inner product. 

It is POSSIBLE TO EXPRESS THE STATE \n) DIRECTLY AS A FUNCTION OF THE GROUND STATE |0) 
USING THE nTH POWER OF THE CREATION OPERATOR AS FOLLOWS: 



|n>: 



■(at) 



t-in- 



/n 



|0> , 



(94) 



which can be obtained by iterations. 

One can also apply this method to get the eigenfunctions in the configuration 

SPACE. To achieve THIS, WE START WITH THE GROUND STATE 



a|0) = 



(95) 



In THE ^-REPRESENTATION, WE HAVE 



ai/)o(a;) : 



x + 



■<po{x) = 



(96) 



Recalling the form of the momentum operator in this representation, we can ob- 
tain A differential equation for the wavefunction of the ground state. Moreover, 

INTRODUCING THE OSCILLATOR LENGTH XQ - 



1 

V'o = . 



WE GET 



(97) 



This equation c.vn be readily solved and the norm.'Vlization to unity of the full line 
integral of the squared modulus of the solution leads to the physical wavefunction 
of the ho ground state 



i'o{x) - 



1 



\f\f^xl 



e 2 ^ xo 



(98) 
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The rest of the eigenfunctions describing the HO excited states, can be obtained 

BY EMPLOYING ITERATIVELY THE CREATION OPERATOR. ThE PROCEDURE IS THE FOLLOWING 
By MATHEMATICAL INDUCTION, ONE CAN SHOW THAT 

5.5 HO Spectrum Obtained from Wilson-Sommerfeld Quan- 
tization Condition 

In the classical phase space, the equation H = E for the HARMONIC oscillator when 

DIVIDED BY E, I.E., 

9 9 9 

+ 5 — = 1 (102) 

2mE 2E ^ ' 

TURNS INTO THE EQUATION FOR AN ELLIPSE 

4+5=1' (103) 

WHERE a = —\— AND h = \j2mE. Therefore, applying the Bohr-Sommerfeld quan- 

TIZATION RULE FOR THIS CASE 

J = (f pxdx = irab = = 2iThn , (104) 

ONE OBTAINS IMMEDIATELY THE SPECTRUM 

En = nhjJo , (105) 
WHICH IS THE QUANTUM HO SPECTRUM UP TO THE ZERO-POINT ENERGY. 



6 The 3D Coulomb Well: The Stationary States 
of the Hydrogen Atom 

The case of the Hydrogen atom corresponds in wave mechanics to an effective 
potential well that is the sum of the coulomb well and the quantum centrifugal 
barrier as shown in fig. (q. this result comes out from the technique of the 
separation of variables that is to be considered for any differential equation in 
more than one dimension. a very good introduction to this technique can be found 

IN THE TEXTBOOK OF ArFKEN AND WeBER [13] . In GENERAL, FOR d VARIABLES THERE ARE 
d — 1 SEPARATION CONSTANTS. 

In SPHERICAL COORDINATES, THE SCHRODINGER EQUATION (V^ + V(r))'4j{r) = Elp{r) 
READS 

ar \ or J sin oti \ aO J sin acp « 

THAT CAN BE ALSO WRITTEN IN THE FORM 

8 f n8ip\ 8 ( d'tb\ d'^ib 2mr'^sm^9 / \ , ^ 

^'^'a-r ['■I) ^^'^'ae [''-'^) ( w = ° ' (^"^^ 

This equation is a partial differential equation for the electron wavefunction 
ip{r,9,(j>) 'within' the atomic hydrogen. Together with the various conditions that 
the wavefunction ip{r, 9, 4>) should fulfill [for example, ■0(r, 9, 4>) SHOULD have a unique 

VALUE AT ANY SPATIAL POINT (r,9,4>)], THIS EQUATION SPECIFIES IN A COMPLETE MANNER THE 
STATIONARY BEHAVIOR OF THE HYDROGEN ELECTRON. 
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n=6,l=2 




n=1,l=0 



-10 







10 
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Figure 7: The effective potential well in the case of the hydrogen atom consisting 
of the electrostatic potential plus a quantized centrifugal barrier (see text). 
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6.1 The Separation of Variables in Spherical Coordinates 



The real usefulness of writing the hydrogen Schrodinger equation in spherical 
coordinates consists in the easy way of achieving the separation procedure in three 
independent one-dimensional equations. the separation procedure is to seek the 
solutions for which the wavefunction 1p{r, d, (j)) has the form of a product of three 
functions, each of one of the three spherical variables, namely r{r) depending only 
on r, 6(5) depending only on 9, and '!'((/)) that depends only on (j). this is quite similar 
to the separation of the laplace equation. thus 

t/'(r, 6», 4)) = fl(r)e(6») $((/)) . (108) 

The R{r) function describes the differential variation of the electron WAVE- 
FUNCTION V ALONG THE VECTOR RADIUS COMING OUT FROM THE NUCLEUS, WITH 9 AND </) 
ASSUMED TO BE CONSTANT. THE DIFFERENTIAL VARIATION OF l/) WITH THE POLAR ANGLE 9 
ALONG A MERIDIAN OF AN ARBITRARY SPHERE CENTERED IN THE NUCLEUS IS DESCRIBED ONLY 
BY THE FUNCTION Q(9) FOR CONSTANT r AND cj). FINALLY, THE FUNCTION $((/)) DESCRIBES 
HOW l/j VARIES WITH THE AZIMUTHAL ANGLE (j) ALONG A PARALLEL OF AN ARBITRARY SPHERE 
CENTERED AT THE NUCLEUS, UNDER THE CONDITIONS THAT r AND 9 ARE KEPT CONSTANT. 

Using i/i = ii0<I>, one can see that 



dr dr 



09 d9 



(109) 



Then, one can obtain the following equations for the three factoring func- 
tions: ^ 

^+mf^ = 0, (110) 



1 

sin ( 



1 d 

f.2 



d 

odR\ 



^de 

d9 



l{l + 1) - 



= 



dr- r 



2m / 



R = 



(111) 



(112) 



h? \47reo'" 

Each of these equations is an ordinary differential equation for a function of 
a single variable. In this way, the Schrodinger equation for the hydrogen elec- 
tron, WHICH INITIALLY WAS A PARTIAL DIFFERENTIAL EQUATION FOR A FUNCTION 1p OF THREE 
VARIABLES, GETS A SIMPLE FORM OF THREE ID ORDINARY DIFFERENTIAL EQUATIONS FOR UN- 
KNOWN FUNCTIONS OF ONE VARIABLE. THE REASON WHY THE SEPARATION CONSTANTS HAVE 
BEEN CHOSEN AS —mf AND —1(1 + 1) WILL BECOME CLEAR IN THE FOLLOWING SUBSECTIONS. 



6.2 The Angular Separation Constants as Quantum Num- 
bers 

6.2.1 The Azimuthal Solution and the Magnetic Quantum Number 

The Eq. lIllOl l is easily solved leading to the following solution 

*(0) =Ar<i,e™'''' , (113) 

where A/# is the integration constant that will be used as a normalization constant 
for the azimuthal part. One of the conditions that any wavefunctions should ful- 
fill IS to have a unique value for any point in space. This applies to <!' as a component 

OF the full WAVEFUNCTION Ip. OnE SHOULD NOTICE THAT <(> AND (j> + 2t7 MUST BE IDENTICAL 
IN THE SAME MERIDIONAL PLANE. THEREFORE, ONE SHOULD HAVE 3>(<^) = <I>(<^ -|- 27r) , I.E., 
J\f^girni4i — J^^^zmi{<l>+2TT) ^ THIS CAN BE FULFILLED ONLY IF m; IS ZERO OR A POSITIVE OR 
NEGATIVE INTEGER (±1, ±2, ±3, ...). THE NUMBER m; IS KNOWN AS THE MAGNETIC QUANTUM 
NUMBER OF THE ATOMIC ELECTRON AND IS RELATED TO THE DIRECTION OF THE PROJECTION OF 
THE ORBITAL MOMENTUM Lz- It COMES INTO PLAY WHENEVER THE EFFECTS OF AXIAL MAG- 
NETIC FIELDS ON THE ELECTRON MAY SHOW UP. THERE IS ALSO A DEEP CONNECTION BETWEEN 
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mi and the orbital quantum number /, which in turn determines the modulus of the 
orbital momentum of the electron. 

The solution for $ should also fulfill the normalization condition when inte- 
grating OVER A full period OF THE AZIMUTHAL ANGLE, 



AND SUBSTITUTING 3>, ONE GETS 



27r 

I <!> |2 d0 = 1 (114) 





Afidil) = 1 . (115) 



It follows that A/* = l/v27r, and therefore the normalized $ IS 



$(0) = ^Le""!'?' . (116) 



6.2.2 The Polar Solution and the Orbital Quantum Number 

The solution of the Q{9) equation is more complicated since it contains two sepa- 
ration CONSTANTS WHICH CAN BE PROVED TO BE INTEGER NUMBERS. THINGS GET EASIER IF 
ONE REMINDS THAT THE SAME EQ. (Illlll OCCURS ALSO WHEN THE HeLMHOLTZ EQUATION FOR 
THE SPATIAL AMPLITUDE PROFILES OF THE ELECTROMAGNETIC NORMAL MODES IS SEPARATED IN 
SPHERICAL COORDINATES. FROM THIS CASE WE ACTUALLY KNOW THAT THIS EQUATION IS THE 
ASSOCIATED LeGENDRE EQUATION FOR WHICH THE POLYNOMIAL SOLUTIONS ARE THE ASSOCI- 
ATED Legendre polynomials 

The function Q{9) are a normalized form, 0(6) = A/eP™', of the associated Legen- 
dre polynomials 

r^iQ\ I '21 + I (I- ^mi, 

For the purposes here, the most important property of these functions is that 
they exist only when the constant I IS an integer number greater or at least equal 

TO \mi\, WHICH IS the absolute VALUE OF m; . THIS CONDITION CAN BE WRITTEN IN THE 
FORM OF THE SET OF VALUES AVAILABLE FOR m,; 

mi = 0,±1,±2,...,±« . (119) 

The CONDITION THAT I SHOULD BE A POSITIVE INTEGER CAN BE SEEN FROM THE FACT THAT 
FOR NONINTEGER VALUES, THE SOLUTION OF EQ. (Illll) DIVERGES FOR COS 6 = ±1, WHILE FOR 
PHYSICAL REASONS WE REQUIRE FINITE SOLUTIONS IN THESE LIMITS. THE OTHER CONDITION 
I > \mi\ CAN BE OBTAINED FROM EXAMINING EQ. (111711 . WHERE ONE CAN SEE A DERIVATIVE 
OF ORDER mi +1 APPLIED TO A POLYNOMIAL OF ORDER 21. THUS m; +1 CANNOT BE GREATER 
THAN 21. On the other hand, the derivative of NEGATIVE ORDER ARE NOT DEFINED AND 
THAT OF ZERO ORDER IS INTERPRETED AS THE UNIT OPERATOR. THIS LEADS TO mi > —I. 

The INTERPRETATION OF THE ORBITAL NUMBER I DOES HAVE SOME DIFFICULTIES. LET US 
EXAMINE THE EQUATION CORRESPONDING TO THE RADIAL WAVEFUNCTION R{r) . THIS EQUATION 
RULES ONLY THE RADIAL MOTION OF THE ELECTRON, I.E., WITH THE RELATIVE DISTANCE WITH 
RESPECT TO THE NUCLEUS ALONG SOME GUIDING ELLIPSES. HOWEVER, THE TOTAL ENERGY OF 
THE ELECTRON E IS ALSO PRESENT. THIS ENERGY INCLUDES THE KINETIC ELECTRON ENERGY IN 
ITS ORBITAL MOTION THAT IS NOT RELATED TO THE RADIAL MOTION. THIS CONTRADICTION CAN 
BE ELIMINATED USING THE FOLLOWING ARGUMENT. THE KINETIC ENERGY T HAS TWO PARTS: 
A PURE RADIAL ONE T^adial AND TgrbUal, WHICH IS DUE TO THE CLOSED ORBITAL MOTION. 
The POTENTIAL ENERGY V OF THE ELECTRON IS THE ATTRACTIVE ELECTROSTATIC ENERGY. 

Therefore, the electron total energy is 

E = Tradial + ^orbital — ■ (120) 

47reo^ 
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Substituting this expression of B in Eq. II112I I we get after some regrouping of the 

TERMS 



1 d 




2m 


y.2 ^j, 







_ hH{i + 1) 

' radial ^ ^ orbital „ 9 



R = Q . (121) 



If the last two terms in parentheses compensate each other, we get a DIFFER- 
ENTIAL equation for THE PURE RADIAL MOTION. THUS, WE IMPOSE THE CONDITION 

J-orbital — ' 5 • i-l-^^j 

However, the orbital kinetic energy of the electron is TorbHal = ^''^'"trbital '^^^ 
SINCE THE ORBITAL MOMENTUM OF THE ELECTRON IS L = raVorbitalT , WE CAN EXPRESS THE 
ORBITAL KINETIC ENERGY IN THE FORM 



Torbital — T 7; ■ (123) 



Therefore, we have 



AND consequently 



L2 _ h^l(l + l) 



2mr-^ 2mr2 



(124) 



L = ^/l{l + l)h . (125) 
The interpretation of this result is that since the orbital quantum number I 

IS constrained to take the values / = 0, 1, 2, (n — 1), THE ELECTRON CAN ONLY HAVE 
ORBITAL MOMENTA SPECIFIED BY MEANS OF EQ. I|125|l . AS IN THE CASE OF THE TOTAL ENERGY 
E, THE ANGULAR MOMENTUM IS CONSERVED AND GETS QUANTIZED. ITS NATURAL UNIT IN 
QUANTUM MECHANICS IS H = h/2-K = 1.054 X 10~^* J.S. 

In the MACROSCOPIC PLANETARY MOTION (PUTTING ASIDE THE MANY-BODY FEATURES), 
THE ORBITAL QUANTUM NUMBER IS SO LARGE THAT ANY DIRECT EXPERIMENTAL DETECTION OF 
THE QUANTUM ORBITAL MOMENTUM IS IMPOSSIBLE. FOR EXAMPLE, AN ELECTRON WITH I = 2 
HAS AN ANGULAR MOMENTUM L = 2.6 X lO"^"* J.S., WHEREAS THE TERRESTRIAL ANGULAR 
MOMENTUM IS 2.7 X 10*" J.S.! 

A COMMON NOTATION FOR THE ANGULAR MOMENTUM STATES IS BY MEANS OF THE LETTER 
S FOR I = 0, p FOR I = 1, d FOR I = 2, AND SO ON. THIS ALPHABETIC CODE COMES FROM THE 
EMPIRICAL SPECTROSCOPIC CLASSIFICATION IN TERMS OF THE SO-CALLED SERIES, WHICH WAS 
IN USE BEFORE THE ADVENT OF WAVE MECHANICS. 

On the other hand, for the interpretation of the MAGNETIC QUANTUM NUMBER, 
WE MUST TAKE INTO ACCOUNT THAT THE ORBITAL MOMENTUM IS A VECTOR OPERATOR AND 
THEREFORE ONE HAS TO SPECIFY ITS DIRECTION, SENSE, AND MODULUS. L, BEING A VECTOR, 
IS PERPENDICULAR ON THE PLANE OF ROTATION. THE GEOMETRIC RULES OF THE VECTORIAL 
PRODUCTS STILL HOLD, IN PARTICULAR THE RULE OF THE RIGHT HAND: ITS DIRECTION AND 
SENSE ARE GIVEN BY THE RIGHT THUMB WHENEVER THE OTHER FOUR FINGERS POINT AT THE 
DIRECTION OF ROTATION. 



6.2.3 The Space Quantization 

We have ALREADY SEEN THE SPATIAL QUANTIZATION OF THE BOHR-SOMMERFELD ELECTRON 
TRAJECTORIES. BUT WHAT SIGNIFICANCE CAN BE ASSOCIATED TO A DIRECTION AND SENSE IN 
THE LIMITED VOLUME OF THE ATOMIC HYDROGEN IN SCHRODINGER WAVE MECHANICS ? THE 
ANSWER MAY BE QUICK IF WE THINK THAT THE ROTATING ELECTRON IS NOTHING BUT A ONE- 
ELECTRON LOOP CURRENT THAT CONSIDERED AS A MAGNETIC DIPOLE HAS A CORRESPONDING 
MAGNETIC FIELD. CONSEQUENTLY, AN ATOMIC ELECTRON WILL ALWAYS INTERACT WITH AN 
APPLIED MAGNETIC FIELD H. THE MAGNETIC QUANTUM NUMBER mj SPECIFIES THE SPATIAL 
DIRECTION OF L, WHICH IS DETERMINED BY THE COMPONENT OF L ALONG THE DIRECTION OF 
THE EXTERNAL MAGNETIC FIELD. THIS EFFECT IS COMMONLY KNOWN AS THE QUANTIZATION 
OF THE SPACE IN A MAGNETIC FIELD. 
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If we choose the direction of the magnetic field as the z axis, the component 

OF L ALONG this DIRECTION IS 

= mih . (126) 
The possible values of m; for a given value of /, go from +1 to — /, passing through 

ZERO, so that there ARE 21 + 1 POSSIBLE ORIENTATIONS OF THE ANGULAR MOMENTUM L IN 
A MAGNETIC FIELD. WHEN I = 0, CAN BE ONLY ZERO; WHENZ = 1, CAN BE H, 0, OR 

— h; WHEN I = 2, Lz TAKES ONLY ONE OF THE VALUES 2h, H, 0, —H, OR —2h, AND SO FORTH. 
It is worth MENTIONING THAT L CANNOT BE PUT EXACTLY PARALL EL OR ANTI-PARALLEL TO 
H, BECAUSE Lz IS ALWAYS SMALLER THAN THE MODULUS + l)h OF THE TOTAL ORBITAL 

MOMENTUM. 

One SHOULD consider the ATOM/eLECTRON CHARACTERIZED BY A GIVEN m; AS HAVING 
THE ORIENTATION OF ITS ANGULAR MOMENTUM L DETERMINED RELATIVE TO THE EXTERNAL 
APPLIED MAGNETIC FIELD. 

In THE ABSENCE OF THE EXTERNAL MAGNETIC FIELD, THE DIRECTION OF THE Z AXIS IS 
FULLY ARBITRARY. THEREFORE, THE COMPONENT OF L IN ANY ARBITRARY CHOSEN DIRECTION 
IS mih; THE EXTERN.'VL MAGNETIC FIELD OFFERS A PREFERRED REFERENCE DIRECTION FROM 
THE EXPERIMENTAL VIEWPOINT. 

Why ONLY THE COMPONENT IS QUANTIZED ? THE ANSWER IS REL;\TED TO THE FACT 
THAT L C.'VNNOT BE PUT ALONG A DIRECTION IN AN ARBITRARY WAY. THERE IS A SPECIAL PRE- 
CESSIONAL MOTION IN WHICH ITS 'VECTORIAL ARROW' MOVES ALWAYS ,'\LONG A CONE CENTERED 
ON THE QUANTIZATION AXIS SUCH THAT ITS PROJECTION Lz IS mih. THE REASON WHY THIS 
QUANTUM PRECESSION OCCURS IS DIFFERENT FROM THE MAC;R0SC0PIC PLANETARY MOTION AS 
IT IS DUE TO THE UNCERTAINTY PRINCIPLE. If L WOULD BE FIXED IN SPACE, IN SUCH A WAY 
THAT Lx, Ly AND Lz WOULD HAVE WELL-DEFINED VALUES, THE ELECTRON WOULD HAVE TO 
BE CONFINED TO A WELL-DEFINED PLANE. FOR EXAMPLE, IF L WOULD BE FIXED ALONG THE Z 
DIRECTION, THE ELECTRON TENDS TO MAINTAIN ITSELF IN THE PLANE xy. 

This can only occur in the case in which the component pz of the electron 

MOMENTUM IS 'INFINITELY' UNCERTAIN. THIS IS HOWEVER IMPOSSIBLE IF THE ELECTRON IS PART 
OF THE HYDROGEN ATOM. BUT SINCE IN REALITY JUST THE COMPONENT Lz OF L TOGETHER 
WITH HAVE WELL DEFINED VALUES AND |L| > \Lz\, THE ELECTRON IS NOT CONSTRAINED 
TO A SINGLE PLANE. If THIS WOULD BE THE CASE, AN UNCERTAINTY WOULD EXIST IN THE 
COORDINATE Z OF THE ELECTRON. THE DIRECTION OF L CHANGES CONTINUOUSLY SO THAT THE 
MEAN VALUES OF Lx AND Ly ARE ZERO, ALTHOUGH Lz KEEPS ALL THE TIME ITS VALUE miH. It 
IS HERE WHERE HeISENBERG'S UNCERTAINTY PRINCIPLE HELPS TO MAKE A CLEAR DIFFERENCE 
BETWEEN ATOMIC WAVE MOTION AND BOHR-SOMMERFELD QUANTIZED ELLIPSES. 

6.3 Polar and Azimuthal Solutions Set Together 

The SOLUTIONS of the azimuthal and POLAR PARTS C;AN be unified WITHIN SPHERICAL 
HARMONICS FUNCTIONS THAT DEPEND ON BOTH (p AND 9. THIS SIMPLIFIES THE ALGEBRAIC 
MANIPULATIONS OF THE FULL WAVE FUNCTIONS tp{r,d,4>). SPHERICAL HARMONICS ARE GIVEN 
BY 

vp {e, 0) = (-1)-' .^±ll^^pp {cose)e'"^"f . (127) 
^ 4-K (1 + mi)\ 

The factor (— 1)*"' does not produce any problem because the Schrodinger equation 
IS linear and homogeneous. This factor is added for the sake of convenience in 

ANGULAR MOMENTUM STUDIES. It IS KNOWN AS THE CONDON-SHORTLEY PHASE FACTOR AND 
ITS EFFECT IS TO INTRODUCE THE ALTERNATED SEQUENCE OF THE SIGNS ± FOR THE SPHERICAL 
HARMONICS OF A GIVEN /. 

6.4 The Radial Solution and the Principal Quantum Num- 
ber 

There is no energy parameter in the angular equations and that is why the angular 

MOTION does not MAKE ANY CONTRIBUTION TO THE HYDROGEN SPECTRUM. It IS THE RADIAL 
MOTION THAT DETERMINES THE ENERGY EIGENVALUES. THE SOLUTION FOR THE RADIAL PART 
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R{r) OF THE WAVE FUNCTION 1p OF THE HYDROGEN ATOM IS SOMEWHAT MORE COMPLICATED 
ALTHOUGH THE PRESENCE OF TWO SEPARATION CONSTANTS, E AND I, POINT TO SOME ASSO- 
CIATED ORTHOGONAL POLYNOMIALS. In THE RADIAL MOTION OF THE HYDROCJEN ELECTRON 
SIGNIFICANT DIFFERENCES WITH RESPECT TO THE ELECTROSTATIC LAPLACE EQUATION DO OC- 
CUR. The flxal result is expressed analytically in terms of the associated Laguerre 

POLYNOMLALS (SCJIIRODINGER 1926). TlIE RADIAL EQUATION CAN BE SOLVED EXACTLY ONLY 
WHEN E IS POSITIVE OR FOR ONE OF THE FOLLOWING NECJATIVE VALUES En (iN WHICH CASES, 
THE ELECTRON IS IN A BOUND STATIONARY STATE WITHIN ATOMIC HYDROGEN) 



-Ry 



WHERE Rv = "i'^'^^o = 13.606 eV IS THE RYDBERG ATOMIC ENERGY UNIT CONNECTED WITH 

THE SPECTROSCOPIC RyDBERG CONSTANT Roo THROUGH RY = hcRoo , WHEREAS nIS A POSITIVE 
INTEGER NUMBER CALLED THE PRINCIPAL QUANTUM NUMBER. It GIVES THE QUANTIZATION OF 
THE ELECTRON ENERGY IN THE HYDROGEN ATOM. THIS DISCRETE ATOMIC SPECTRUM HAS BEEN 
FIRST OBTAINED IN 1913 BY BOHR USING SEMI-EMPIRICAL QUANTIZATION METHODS AND NEXT 
BY PAULI AND SCHRODINGER ALMOST SIMULTANEOUSLY IN 1926. 

Another condition that should be satisfied in order to solve the radial equation 

IS that n HAVE TO BE STRICTLY BIGGER THAN I. ItS LOWEST VALUE IS / -f 1 FOR A GIVEN I. 

Vice versa, the condition on / is 

« = 0,l,2,...,(n-l) , (129) 

FOR GIVEN n. 

The radial equation can be written in the form 



„2 



d'^R dR 

— 5- +2r 

dr-^ dr 



2mE , 2me2 



- l{l + 1) 



R = Q . (130) 



Dividing by and usinc; the substitution x(''') = i'R to eliminate the first derivative 

dR 

dr 



— ONE GETS THE STAND.ARD FORM OF THE RADIAL SCHRODINGER EQUATION DISPLAYING THE 



effective POTENTIAL U [r) = — COIlst/r + i(/ + l)/r-^ (.actually, ELECTROSTATIC POTENTIAL 
PLUS QUANTIZED CENTRIFUGAL BARRIER). THESE ARE NECESSARY MATHEMATICAL STEPS IN 
ORDER TO DISCUSS A NEW BOUNDARY CONDITION, SINCE THE SPECTRUM IS OBTAINED BY MEANS 
OF THE R EQUATION. THE DIFFERENCE BETWEEN A RADIAL SCHRODINGER EQUATION AND 
A FULL-LINE ONE IS THAT A SUPPLEMENTARY BOUNDARY CONDITION SHOULD BE IMPOSED AT 

THE ORIGIN (r = 0). The Coulomb potential belongs to a class of potentials that 

ARE called WEAK SINGULAR FOR WHICH lim,.^o U{r)r'^ = 0. IN THESE CASES, ONE TRIES 
SOLUTIONS OF THE TYPE X r", IMPLYING I/(l/ — 1) = l{l + 1), SO THAT THE SOLUTIONS ARE 
Ul =1 + 1 AND 1/2 = —I, JUST AS IN ELECTROSTATICS. THE NEGATIVE SOLUTION IS ELIMINATED 
FOR I 7^ BECAUSE IT LEADS TO A DIVERGENT NORMALIZATION CONSTANT, NOR DID IT RESPECT 
THE NORMALIZATION AT THE DELTA FUNCTION FOR THE CONTINUOUS PART OF THE SPECTRUM. 
On THE OTHER HAND, THE PARTICULAR CASE J/2 = IS ELIMINATED BECAUSE THE MEAN KINETIC 

energy is not finite. the final conclusion is that x(0) = for any i. 

Going back to the analysis of the radial equation for R, the first thing to do 
IS to write it in nondimensional variables. This is performed by noticing tilat the 
only space and time scales th.at one can form on combininc; the three FL'NDAMENTAL 

constants entering this PROBLEM, N.AMELY e^, irie AND ft ARE THE BOHR ATOMIC RADIUS 

as = f? /me^ = 0.529 • 10^* cm. and the atomic time ts = f? jvae^ = 0.242 • 10~^^ SEC, 

USUALLY KNOWN AS ATOMIC UNITS. EMPLOYING THESE UNITS, ONE GETS 



fR. 2dR 
dr^ r dr 



R = 0, (131) 



WHERE WE ARE ESPECIALLY INTERESTED IN THE DISCRETE PART OF THE SPECTRUM (E <0). 

The NOTATIONS n = l/\/—E and p = 2r/n leads us to 



d'^R 2 dR 
\ . 

dp^ p dp 



n _1 _l{l + l) 
p 4 p2 



R = . (132) 
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For p — > oo, this equation reduces to = -j, having solutions R <x e^i'l'^ . Be- 
cause OF THE normalization CONDITION ONLY THE DECAYING EXPONENTIAL IS ACCEPTABLE. 

On the other hand, the asymptotic at zero, as we already commented ON, should 
BE _R oc p'. Therefore, we can write _R as a product of three radial functions 
_R = p'e~''/^F(p), OF which the first two give the asymptotic behaviors, whereas the 

THIRD IS THE RADIAL FUNCTION IN THE INTERMEDIATE REGION. THE LATTER FUNCTION IS OF 
MOST INTEREST BECAUSE ITS FEATURES DETERMINE THE ENERGY SPECTRUM. THE EQUATION 
FOR F IS 

(p- F dF 
p—^ + {2l + 2-p)— + {n-l-l)F = 0. (133) 
dp'' dp 

This is a particular case of confluent hypergeometric equation. It can be identi- 
fied AS THE EQUATION FOR THE ASSOCIATED LAGUERRE POLYNOMIALS L^^^{p). THUS, THE 
NORMALIZED FORM OF R IS 



RM = -4J o^ ' 'li e-^/V^^'+VW . (134) 

n/= y 2n[{n + l)]]-^ 



WHERE THE FOLLOWING LaGUERRE NORMALIZATION CONDITION HAS BEEN USED 

-V-ML--(P)]^.P=|[^^^. (135) 
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6.5 Final Formulas for the Hydrogen Atom Stationary 
States 

We HAVE NOW THE SOLUTIONS OF ALL THE EQUATIONS DEPENDING ON A SINGLE VARIABLE AND 
THEREFORE WE CAN BUILD THE STATIONARY WAVE FUNCTIONS FOR ANY ELECTRONIC STATE OF 
THE HYDROGEN ATOM. THEY HAVE THE FOLLOWING ANALYTIC FORM 

4^(r,e,4,)=MHianrye-°'"'-^^L^J+\anr)P;"'icose)e'"''1' , (136) 



2 21 + 1 (l-mi)\ (n-l-iy. 
y 47r (l + mi)\ [{n + l)\\-^ 
Using the spherical harmonics, the solution is written as follows 



If the dynamical factor is included, we get 



V.(r,e,</,;t) = -A^(!!_J_J^(^^ . (139) 

The latter formulas may be considered as the final result for the Schrodinger 

SOLUTION OF the HYDROGEN ATOM FOR ANY STATIONARY ELECTRON STATE. INDEED, ONE CAN 
SEE EXPLICITLY BOTH THE ASYMPTOTIC DEPENDENCE AND THE TWO ORTHOGONAL AND COM- 
PLETE SETS OF FUNCTIONS, I.E., THE ASSOCIATED LaGUERRE POLYNOMIALS AND THE SPHERICAL 
HARMONICS THAT CORRESPOND TO THIS PARTICULAR CASE OF LINEAR PARTIAL SECOND-ORDER 
DIFFERENTIAL EQUATION. FOR THE ALGEBRAIC APPROACH TO THE HYDROGEN ATOM PROBLEM 
WE RECOMMEND THE PAPER OF KiRCHBERG AND COLLABORATORS [16] AND THE REFERENCES 

THEREIN. Finally, the stationary hydrogen eigenfunctions are characterized by a 

HIGH DEGENERACY SINCE THEY DEPEND ON THREE QUANTUM NUMBERS WHEREAS THE ENERGY 
spectrum IS ONLY n-DEPENDENT. THE DEGREE OF DEGENERACY IS EASILY CALCULATED IF ONE 
NOTICE THAT THERE ARE 21 + 1 VALUES OF m; FOR A GIVEN I WHICH IN TURN TAKES VALUES 
FROM TO n — 1 FOR A GIVEN n. THUS, THERE ARE J2"^o C^^ + 1) = n'^ STATES OF GIVEN 
ENERGY. Not only THE PRESENCE OF MANY DEGREES OF FREEDOM IS THE CAUSE OF THE 
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STRONG DEGENERACY. In THE CASE OF THE HYDROGEN ATOM THE EXISTENCE OF THE CON- 
SERVED Runge-Lenz vector (commuting operator)K = {2me^)~^[L, X P — P X L] + r/r 

INTRODUCES MORE SYMMETRY INTO THE PROBLEM AND ENHANCES THE DEGENERACY. On THE 
OTHER HAND, THE ONE-DIMENSIONAL QUANTUM WAVEFUNCTIONS ARE NOT DEGENERATE BE- 
CAUSE THEY ARE CHARACTERIZED BY A SINGLE DISCRETE INDEX. THE GENERAL PROBLEM OF 
DEGENERACIES IS NICELY PRESENTED IN A PAPER BY ShEA AND ArAVIND [TTl . 

6.6 Electronic Probability Density 

In THE Bohr model of the hydrogen atom, the electron rotates around the nucleus 

ON CIRCULAR OR ELLIPTIC TRAJECTORIES. It IS POSSIBLE TO THINK OF APPROPRIATE EXPERI- 
MENTS allowing TO "see" that THE ELECTRON MOVES WITHIN EXPERIMENTAL ERRORS AT THE 
PREDICTED RADII Vn = n?ao IN THE EQUATORIAL PLANE 6 = 90°, WHEREAS THE AZIMUTHAL 
ANGLE MAY VARY ACCORDING TO THE SPECIFIC EXPERIMENTAL CONDITIONS. 

It is in this case that the more general WAVE MECHANICS CHANGES THE CONCLUSIONS 
OF THE Bohr model in at least two important ASPECTS: 

• First, one cannot speak about exact values of r, d, 4> (and therefore of planetary 
trajectories), but only of relative probabilities to find the electron within a given 
infinitesimal region of space. This feature is a consequence of the wave nature of 
the electron. 

• Secondly, the electron does not move around the nucleus in the classical con- 
ventional WAY because the PROBABILITY DENSITY \tl>\^ DOES NOT DEPEND ON TIME BUT 
CAN VARY SUBSTANTIALLY AS A FUNCTION OF THE RELATIVE POSITION OF THE INFINITESIMAL 
REGION. 

The hydrogen electron wave function ^ is tp = RO^, where R = Rni{r) describes 

THE WAY l/j changes WITH r WHEN THE PRINCIPAL AND ORBITAL QUANTUM NUMBERS HAVE THE 
VALUES n AND /, RESPECTIVELY. = Qlm, (0) DESCRIBES IN TURN HOW 1p VARIES WITH 9 WHEN 
THE ORBITAL AND MAGNETIC QUANTUM NUMBERS HAVE THE VALUES I AND m; , RESPECTIVELY. 

Finally, 3> = <E>m( (</>) gives the change of ^ with </> when the magnetic quantum 

NUMBER HAS THE VALUE m; . THE PROBABILITY DENSITY \ 1p CAN BE WRITTEN 

|^|2 = |R|2|e|2|<l>|2 . (140) 

Notice that the probability density which measures the possibility to find 
the electron at a given azimuthal angle </>, is a constant (does not depend on 4>) ■ 
Therefore, the electronic probability density is symmetric with respect to the z 

AXIS AND INDEPENDENT ON THE MAGNETIC SUBSTATES (aT LEAST UNTIL AN EXTERNAL MAGNETIC 

field is applied) . Consequently, the electron has an equal probability to be found 

IN any AZIMUTHAL DIRECTION. THE RADIAL PART R OF THE WAVE FUNCTION, CONTRARY TO 
NOT ONLY VARIES WITH r, BUT IT DOES IT DIFFERENTLY FOR ANY DIFFERENT COMBINATION 
OF QUANTUM NUMBERS n AND I. FIGURE ^ SHOWS PLOTS OF J? AS A FUNCTION OF r FOR 
THE STATES Is, 2s, AND 2p. R IS MAXIMUM AT THE CENTER OF THE NUCLEUS (r = 0) FOR 
ALL THE S STATES, WHEREAS IT IS ZERO AT r = FOR ALL THE STATES OF NONZERO ANGULAR 
MOMENTUM. 

The ELECTRON PROBABILITY DENSITY AT THE POINT r,9,4> IS PROPORTIONAL TO IV'P, BUT 
THE REAL PROBABILITY IN THE INFINITESIMAL VOLUME ELEMENT dV IS | l/" P dV . In SPHERICAL 
COORDINATES dV = r'^ sin 9 drdBdlj) AND SINCE AND # ARE NORMALIZED FUNCTIONS, THE 
REAL NUMERICAL PROBABILITY P(r)dr TO FIND THE ELECTRON AT A RELATIVE DISTANCE WITH 
RESPECT TO THE NUCLEUS BETWEEN r AND r + dr IS 

P{r)dr = \ R]"^ dr | 1^ sined9 / \ <S> \^ d<j> 
Jo Jo 

= \ R\'^ dr . (141) 
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Figure 9: Probability density of finding the liydrogen electron between r and 
r + dr with respect to the nucleus for the states Is (blue), 2s (red), 2p (green). 
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The function P(r) is displayed in Fig. ([9j for the same states for which the 

RADIAL FUNCTIONS R IS DISPLAYED IN FiG. ((Sjl . In PRINCIPLE, THE CURVES ARE QUITE DIF- 
FERENT. We immediately see that P(r) is not maximal in the nucleus for the states 

S, AS HAPPENS FOR R. INSTEAD, THEIR MAXIMA ARE ENCOUNTERED AT A FINITE DISTANCE 

from the nucleus. The most probable value of r for a Is electron is exactly as, 
THE Bohr radius. However, the mean value of r for a Is electron is l.bag. At 

FIRST sight this MIGHT LOOK STRANGE, BECAUSE THE ENERGY LEVELS ARE THE SAME BOTH 
IN QUANTUM MECHANICS AND IN BOHR'S MODEL. THIS APPARENT OUTMATCHING IS ELIMINATED 
IF ONE TAKES INTO ACCOUNT THAT THE ELECTRON ENERGY DEPENDS ON 1/r AND NOT ON r, 
AND THE MEAN VALUE OF 1/r FOR A Is ELECTRON IS EXACTLY l/aQ. 

The function varies with the polar angle 9 for all the quantum numbers I 

AND nil, UNLESS Z = mj = 0, WHICH ARE THE S STATES. THE PROBABILITY DENSITY |0p 
FOR A S STATE IS A CONSTANT (1/2). THIS MEANS THAT SINCE |<I>P IS ALSO A CONSTANT, 
THE ELECTRONIC PROBABILITY DENSITY HAS THE SAME VALUE FOR A GIVEN VALUE OF r, 

NOT DEPENDING ON THE DIRECTION. In OTHER STATES, THE ELECTRONS PRESENT AN ANGULAR 
BEHAVIOR THAT IN MANY CASES MAY BE QUITE COMPLICATED. BECAUSE \lp\'^ IS INDEPENDENT 
OF A THREE-DIMENSIONAL REPRESENTATION OF \tl>\^ CAN BE OBTAINED BY ROTATING A 
PARTICULAR REPRESENTATION AROUND A VERTICAL AXIS. THIS CAN PROVE VISUALLY THAT 
THE PROBABILITY DENSITIES FOR THE S STATES HAVE SPHERICAL SYMMETRY, WHILE ALL THE 
OTHER STATES DO NOT POSSESS IT. In THIS WAY, ONE CAN GET MORE OR LESS PRONOUNCED 
LOBES OF CHARACTERISTIC FORMS DEPENDING ON THE ATOMIC STATE. THESE LOBES ARE QUITE 
IMPORTANT IN CHEMISTRY FOR SPECIFYING THE ATOMIC INTERACTION IN THE MOLECULAR BULK. 



6.7 Other 3D Coordinate Systems Allowing Separation of 
Variables 

a complete discussion of the 3d coordinate systems allowing the separation of 
variables for the schrodinger equation has been provided by cook and fowler i18| . 
Here, we briefly review these coordinate systems: 
Parabolic. 

The parabolic coordinates given by 

X = cos , 

y = \/???sin0 , (142) 



WHERE 5 e [0, oo), r] e [0, oo), AND (/> S [0, 27r), ARE ANOTHER COORDINATE SYSTEM IN 
WHICH THE SCHRODINGER HYDROGEN EQUATION IS SEPARABLE AS FIRST SHOWN BY SCHRODINGER 
[14] . The final solution in this case is expressed as the product of factors OF ASYMP- 
TOTIC NATURE, AZIMUTHAL HARMONICS, AND TWO SETS OF ASSOCIATE LAGUERRE POLYNOMIALS 
IN THE VARIABLES J AND Tj, RESPECTIVELY. THE ENERGY SPECTRUM (— l/n^) AND THE DEGEN- 
ERACY (n^) OF COURSE DO NOT DEPEND ON THE COORDINATE SYSTEM. THEY ARE USUALLY 
EMPLOYED IN THE STUDY OF THE StARK EFFECT AS FIRST SHOWN BY EPSTEIN 1151 . 

Spheroidal. 

Spheroidal coordinates also can be treated by the separation technique with the 

2 COMPONENT OF ANGULAR MOMENTUM REMAINING AS A CONSTANT OF THE MOTION. ThERE 

are two types of spheroidal coordinates. 

The oblate spheroidal coordinates are given by: 

X = r cosh 5 COS -q cos , 

y = r cosh 5 cos ?7 sin iji , (143) 
z = r sinh § sin rj , 

where ^ e [0, oo), 77 e [-7r/2, 7r/2], AND £ [0, 27r). 
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The prolate spheroidal coordinates are complementary to the oblate ones in 

THE variables ^ AND T]: 

X = r sinh § sin r) cos 4> , 

y = r sinh § sin sin , (144) 
z = r cosh 5 cos -q , 
where 5 e [0, oo), rj 6 [0, tt], and cf> e [0, 2tt). 
Spheroconal. 

The spheroconal system is a quite unfamiliar coordinate separation system of the H- 
atom Schrodinger equation. In this case is retained as a constant of the motion 

BUT 1/2 IS replaced BY ANOTHER SEPARATION OPERATOR B. THE RELATIONSHIP WITH THE 

Cartesian system is given through [18] 



(145) 



One of the separation constants is l{l + 1), the eigenvalue of the separation opera- 
tor L-^. This fact is of considerable help in dealing with the unfamiliar equations 

RESULTING FROM THE SEPARATION OF THE 9 AND A EQUATIONS. THE SEPARATION OPERATOR 
B CAN BE TRANSFORMED INTO A CARTESIAN FORM, 

B = b^Ll + C^LI , (146) 

WHERE Ly AND Lz ARE THE USUAL CARTESIAN COMPONENTS OF THE ANGULAR MOMENTUM 
OPERATOR. Thus, in the spheroconal system, a linear combination of the SQUARES 
OF THE y AND Z COMPONENTS OF THE ANGULAR MOMENTUM EFFECTS THE SEPARATION. THE 
LINEAR COMBINATION COEFFICIENTS ARE SIMPLY THE SQUARES OF THE LIMITS OF THE SPHERO- 
CONAL COORDINATE RANGES. 




In GENERAL, THE STATIONARY STATES IN ANY OF THESE COORDINATE SYSTEMS CAN BE 
WRITTEN AS LINEAR COMBINATION OF DEGENERATE EIGENFUNCTIONS OF THE OTHER SYSTEM 
AND THE GROUND STATE (THE VACUUM) SHOULD BE THE SAME. 



7 The 3D Parabolic Well: The Stationary States 
of the Isotropic Harmonic Oscillator 

We COMMENTED ON THE IMPORTANCE IN PHYSICS OF THE HO AT THE BEGINNING OF OUR 
ANALYSIS OF THE ID QUANTUM HO. If WE WILL CONSIDER A 3D ANALOG, WE WOULD BE LED 

TO STUDY A Taylor expansion of the potential in all three variables retaining the 

TERMS UP TO the SECOND ORDER, WHICH IS A QUADRATIC FORM OF THE MOST GENERAL FORM 

V{x, y, z) = ax^ -|- by^ + cz^ + dxy + exz + fyz . (147) 
There are however many systems with spherical symmetry or for which this 

SYMMETRY IS SUFFICIENTLY EXACT. ThEN, THE POTENTIAL TAKES THE MUCH SIMPLER FORM 

V{x,y,z) = K{x'^ + y^ + z'^) . (148) 

This is equivalent to assuming that the second unmixed partial spatial derivatives 
OF the potential have all the same value, herein denoted BY' K. We can add up 
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that this is a good approximation whenever the values of the mixed second partial 
derivatives are small in comparison to the unmixed ones. when these conditions 
are satisfied and the potential is given by (i148ii . we say that the system is a 3d 
spherically symmetric ho. 
The Hamiltonian in this case is of the form 

/f= — V'+— r■^ (149) 
2m 2 

where the Laplace operator is given in spherical coordinates and r is the spheri- 
cal RADIAL COORDINATE. EQUIVALENTLY, THE PROBLEM CAN BE CONSIDERED IN CARTESIAN 

coordinates but it is a trivial generalization of the one-dimensional case. 
Since the potential is time independent the energy is conserved. In addition, be- 
cause OF the spherical sy'mmetry the orbital momentum is also conserved. Having 
TWO conserved quantities, we associate to each of them a corresponding quantum 
number. as a matter of fact, as we have seen in the case of the hydrogen atom, the 
spherical symmetry leads to three quantum numbers, but the third one, the mag- 
netic number, is related to the 'space quantization' and not to the geometrical 
features of the motion. 

Thus, the eigenfunctions depend effectively only on two quantum numbers. The 
eigenvalue problem of interest is then 

= E^i^K^i . (150) 
The Laplace operator in spherical coordinates reads 

^2 _ 9^ _l_ 2 9 L'^ ^^^^^ 

where the angular operator is the usual spherical one 

19/ a\ 1 

sin f 



sine de \ de J sine^ d<p^_ 
The eigenfunctions of are the spherical harmonics, i.e. 



(152) 



L^Yirr,, (9, <p) = -h^i(i + i)y,„, (e, i^) . (153) 

In order to achieve the separation of the variables and functions, the following 
substitution is proposed 

*n!m, {r, e, = ^^:^Yi^^ (e, ^) . (i54) 

Once this is plugged in the Schrodinger equation, the spatial and the angular parts 
are separated from one another. the equation for the spatial part has the form 

2rnE„i im?uj'^ ^ 1(1 + 1) 



ft2 ?i2 



r^-^— ^ R„,(r)=0. (155) 



Using the oscillator parameters fc^^ = ^"^f "' and A = the previous equa- 
tion IS precisely of the one-dimensional QUANTUM OSCILLATOR FORM BUT IN THE RADIAL 
VARIABLE AND WITH AN ADDITIONAL ANGULAR MOMENTUM BARRIER TERM, 

, - A^r^ - = . (156) 



To SOLVE THIS EQUATION, WE SHALL START WITH ITS ASYMPTOTIC ANALYSIS. EXAMINING FIRST 
THE INFINITE LIMIT r ^ OO, WE NOTICE THAT THE ORBITAL MOMENTUM TERM IS NEGLIGIBLE 
AND THEREFORE IN THIS LIMIT THE ASYMPTOTIC BEHAVIOR IS SIMILAR TO THAT OF THE ONE- 
DIMENSIONAL OSCILLATOR, I.E., A GAUSSIAN TAIL 

Rnl{r) ^ cxp I ~ — j FOR r ^ OO . (157) 
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If now we consider the behavior near the origin, we can see that the dominant 

TERM is that OF THE ORBITAL MOMENTUM, I.E., THE DIFFERENTIAL EQUATION (I156II IN THIS 
limit turns INTO 

Ki - ^-^^Rni = . (158) 
This is a differential equation of the Euler type 

r"y{n) ^ ^n-iy(n-i) (.r) + . . . + ry'{r) + y(r) = 

FOR THE CASE n = 2 WITH THE FIRST DERIVATIVE MISSING. FOR SUCH EQUATIONS THE SOLU- 
TIONS ARE SOUGHT OF THE FORM y = r°' THAT PLUGGED IN THE EQUATION LEAD TO A SIMPLE 
POLYNOMIAL EQUATION IN a, WHOSE TWO INDEPENDENT SOLUTIONS ARE / + 1 AND — /. THUS, 
ONE GETS 



OR 



FOR 







The PREVIOUS arguments lead to proposing the SUBSTITUTION 

Ar-2- 



Rulir) = r'+^ e^p 



The second possible substitution 



exp 



«(r) 



(159) 
(160) 

(161) 



produces the same equation as II160I I. Substituting II160I I in (I158II . the following 
differential equation for u is obtained 



+ 2 (^'-^ - Ar^ u' - 1(21 + 3)A - kli]u = 



By using now the change of variable to = At-2, one gets 



wu" + {l+--w]u'- 



~2 



, 



(162) 



(163) 



IS THE SAME DIMENSIONLESS ENERGY PARAMETER AS IN THE ONE- 



WHERE l^-nl = ^ = ^ 
DIMENSIONAL CASE BUT NOW WITH TWO SUBINDICES. We SEE THAT WE FOUND AGAIN A DIF- 
FERENTIAL EQUATION OF THE CONFLUENT HYPERGEOMETRIC TYPE HAVING THE SOLUTION 



u{r)=A iFi 



+ B r-(2'+i) iFi 



2 



(164) 



The second particular solution cannot be normalized because it DIVERGES STRONGLY 

FOR r 0. Thus one takes B = which leads to 



1+ Kr^l] ,1+ -\\r 



(165) 



By using the same arguments on the asymptotic BEHAVIOR AS IN THE ONE-DIMENSIONAL 

ho case, that is, imposing a regular solution at infinity, leads to the truncation of 
the confluent hypergeometric series, which implies the quantization of the energy. 
The truncation is explicitly 



AND SUBSTITUTING K„; , WE GET THE ENERGY SPECTRUM 

3 



-n , 



E„i = fiLij [2n + l + 



tuj { N+ ■ 



(166) 



(167) 
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One can notice that for the three-dimensional spherically symmetric HO there is 
A zero point energy of ^fkv, three times bigger than in the one-dimensional case. 
The unnormalized eigenfunctions of the three-dimensional harmonic oscillator are 

V'nim,(r-,e,</p) =r'e(~~) ^Fi (^-n,l + Xr^^ Yi^^{e, ^p) . (168) 
Since we are in the case of a radial problem with centrifugal barrier we know 

THAT WE HAVE TO GET AS SOLUTIONS THE ASSOCIATED LAGUERRE POLYNOMIALS. THIS CAN BE 
SEEN BY USING THE CONDITION II166I I IN (I163II WHICH BECOMES 

wu" + (^l + ^ - tvj u' + nu = . (169) 

The latter equation has the form of the associated Laguerre equation wu" + 
{p + 1 ~ w)u' + nu = for p = 1 + 1/2 with the polynomial solutions l',^'^^^ . Thus, 
the normalized solutions can be written 

V'„,„,(r,e,¥>) =Ar„,r'e("~) (Xr^) Yi,„^{0, <fi) , (170) 

where the normalization integral can be calculated as in the id oscillator case 
using the Laguerre generating function 



-zt/{l-t) fP 
p=0 ' 



The final result is 

~ n\{n + iy. 



(172) 



TT 2-(2"+2!+2)(2n + 21 + 1)!. 

If the dynamical phase factor is included, the stationary wavefunctions of the 
3D spherically symmetric oscillator take the following final form 

V'n.i,™, {r, e, <t>, t) = Xni exp (~i (2n + l + ^ j u,t - r' Ll+^^\\r'') Y;^, (0, ^) . (173) 

For the algebraic (factorization) method applied to the radial oscillator we refer 
the reader to detailed studies [191. The degeneracy of the radial oscillator is 

EASIER TO calculate BY COUNTING THE CARTESIAN EIGENSTATES AT A GIVEN ENERGY, WHICH 

GIVES (Af-|-2)!/Af!2! [17]. 



8 Stationary Bound States in the Continuum 

After all these examples, it seems that potential wells are necessary for the 

EXISTENCE OF STATIONARY STATES IN WAVE MECHANICS. HOWEVER, THIS IS NOT SO! AlL 
OF THE QUANTUM BOUND STATES CONSIDERED SO FAR HAVE THE PROPERTY THAT THE TOTAL 
ENERGY OF THE STATE IS LESS THAN THE VALUE OF THE POTENTIAL ENERGY AT INFINITY, 
WHICH IS SIMILAR TO THE BOUND STATES IN CLASSICAL MECHANICS. THE BOUNDNESS OF THE 
QUANTUM SYSTEM IS DUE TO THE LACK OF SUFFICIENT ENERGY TO DISSOCIATE. HOWEVER, IN 
WAVE MECHANICS IT IS POSSIBLE TO HAVE BOUND STATES THAT DO NOT POSSESS THIS PROPERTY, 
AND WHICH THEREFORE HAVE NO CLASSICAL ANALOG. 

Let us CHOOSE the zero of energy so that THE POTENTIAL ENERGY FUNCTION VANISHES 
AT INFINITY. THE USUAL ENERGY SPECTRUM FOR SUCH A POTENTIAL WOULD BE A POSITIVE 
ENERGY CONTINUUM OF UNBOUND STATES, WITH THE BOUND STATES, IF ANY, OCCURRING AT 
DISCRETE NEGATIVE ENERGIES. HOWEVER, StILLINGER AND HeRRICK (1975) [2D], FOLLOWING 
AN EARLIER SUGGESTION BY VON NEUMANN AND WiGNER [21] , HAVE CONSTRUCTED POTENTIALS 
THAT HAVE DISCRETE BOUND STATES EMBEDDED IN THE POSITIVE ENERGY CONTINUUM. BOUND 
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STATES ARE REPRESENTED BY THOSE SOLUTIONS OF THE EQUATION (— + V)"^ = E<S' FOR 
WHICH THE NORMALIZATION INTEGRAL f\'^>\'^CpX IS FINITE. (We ADOPT UNITS SUCH THAT 
h= 1 AND m =1.) 

We can FORMALLY SOLVE FOR THE POTENTIAL, 

V{r;E) = E+-(^^j . (174) 
For THE POTENTIAL TO BE NONSINGULAR, THE NODES OF * MUST BE MATCHED BY ZEROS OF 

V^^I/. The free particle zero-angular-momentum function 4'o(x) = sm{kr)/kr satis- 
fies II174II with energy eigenvalue Eg = ifc^ AND WITH V IDENTICALLY EQUAL TO ZERO, 
BUT IT IS UNACCEPTABLE BECAUSE THE INTEGRAL OF |^0p IS NOT CONVERGENT. HOWEVER, 
BY TAKING 

vl/(x) = *o(x)/(r) , (175) 

AND REQUIRING THAT /(r) GO TO ZERO MORE RAPIDLY THAN r~^/^ AS r -+ OO ONE CAN GET 
A CONVERGENT INTEGRAL FOR |>I'((x)p. SUBSTITUTING II175I I INTO (I174II . WE OBTAIN 

1 f'(r) 1 f"(r) 

V{r;E) = E--k^ + kcot(kr)l^ + -i-^ . (176) 

2 fir) 2 /(r) 

For V TO REMAIN BOUNDED, f'{r)/f{r) MUST VANISH AT THE POLES OF cot(fer); THAT IS, 
AT THE ZEROS OF sin(fcr). THIS CAN BE ACHIEVED IN DIFFERENT WAYS BUT IN ALL KNOWN 
PROCEDURES THE MODULATION FACTOR /(r) HAS THE FORM 

/(r) = [A + s(r)]-i , (177) 

WHERE A IS A POSITIVE CONSTANT, ALTHOUGH StILLINGER AND HeRRICK MENTIONED A WIDER 
CLASS OF POSSIBLE f{r). THEY CHOSE THE MODULATION VARIABLE 

s^hir) = Sk^ / r'[sm{kr')fdr' = -(2krf - 2fcrsin(2fcr) - cos(2fcr) + 1 . (178) 
Jo 2 

The principles guiding the choice of s(r) are: that the integrand must be nonneg- 

ATIVE, so that s(r) will BE A MONOTONIC FUNCTION OF r; AND THAT THE INTEGRAND MUST 
be PROPORTIONAL TO sin(fcr), SO THAT ds(r)/dr will VANISH AT THE ZEROS OF sin(fer). FOR 
Ssh, * DECREASES LIKE AS r — > OO, WHICH ENSURES ITS SQUARE INTEGRABILITY. THE 

POTENTIAL II176I I THEN BECOMES (fOR. Eq) 

64fe4r2[sin(fcr)]4 4k'^ { [sm(kr)]'^ + 2kr sm(2kr)} 
Vshir;Eo) = —— — — — . (179) 

The ENERGY OF THE BOUND STATE PRODUCED BY THE POTENTIAL V{r;Eo) IS Eq = ifc^ 

as for the free particle, i.e., it is independent of a and the modulation factor f{r). 
Therefore, the main idea for getting bound states in the continuum is to build 
isospectral potentials of the free particle and more generally for any type of 

scattering STATE. A MORE CONSISTENT PROCEDURE TO OBTAIN ISOSPECTRAL POTENTIALS IS 
GIVEN BY THE FORMALISM OF SUPERSYMMETRIC QUANTUM MECHANICS [22] WHICH IS BASED 
ON THE DARBOUX TRANSFORMATIONS . FOR THE SUPERSYMMETRIC CASE THE MODULATION 
VARIABLE IS 

sair)= / u^{T')dr' = / [sm{kT')]'dr' = -r sin(2A;r) (180) 

Jo Jo 2 4k 

AND THE ISOSPECTRAL POTENTIAL HAS THE FORM 

2[sin(fer)l"' 2k sin(2kr) 
Vd(r; Eo) = , ' ^ / - . ^ \ ' ■ (181) 

Finally, in the amplitude modulation method of von Neumann and Wigner the mod- 
ulation VARIABLE IS GIVEN BY 

s„„„(r) = (4fcsd)^ = [2fcr - sin(2fcr)]^ (182) 
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AND THE ISOSPECTRAL FREE PARTICLE POTENTIAL IS 




(183) 



Vshir;Eo) ~ Va{r;Eo) ~ Vv„^ir;Eo) 5^ '- = ~Sk^smc(2kr) = 




8k^jo{2kr) 



(184) 



WHERE THE SINC IS THE CARDINAL SINE FUNCTION WHICH IS IDENTICAL TO THE SPHERICAL 
BeSSEL FUNCTION OF THE FIRST KIND JQ- On THE OTHER HAND, THESE POTENTIALS DISPLAY 

different power-law behavior near the origin. 

Moreover, the asymptotic sinc form which is typical in diffraction suggests 
that the existence of bound states in the continuum can be understood by using the 

analogy of WAVE PROPAGATION TO DESCRIBE THE DYNAMICS OF QUANTUM STATES. It SEEMS 
THAT THE MECHANISM WHICH PREVENTS THE BOUND STATE FROM DISPERSING LIKE ORDINARY 
POSITIVE ENERGY STATES IS THE DESTRUCTIVE INTERFERENCE OF THE WAVES REFLECTED FROM 
THE OSCILLATIONS OF V{r; E). ACCORDING TO StILLINGER AND HeRRICK NO OTHER f(r) THAT 
PRODUCES A SINGLE PARTICLE BOUND STATE IN THE CONTINUUM WILL LEAD TO A POTENTIAL 
THAT DECAYS MORE RAPIDLY THAN l|184|l . HOWEVER, THEY PRESENT FURTHER DETAILS IN THEIR 
PAPER WHICH SUGGEST THAT NONSEPARABLE MULTIPARTICLE SYSTEMS, SUCH AS TWO- ELECTRON 
ATOMS, MAY POSSESS BOUNDED STATES IN THE CONTINUUM WITHOUT SUCH A CONTRIVED FORM 
OF POTENTIAL AS (|179p . 

The BOUND (or localized) states in THE CONTINUUM ARE NOT ONLY OF ACADEMIC 
INTEREST. Such a type of electronic stationary QUANTUM STATE HAS BEEN PUT INTO 
EVIDENCE IN 1992 BY CAPASSO AND COLLABORATORS [24] BY INFRARED ABSORPTION MEA- 
SUREMENTS ON A SEMICONDUCTOR SUPERLATTICE GROWN BY MOLECULAR BEAM EPITAXY IN 
SUCH A WAY THAT ONE THICK QUANTUM WELL IS SURROUNDED ON BOTH SIDES BY SEVERAL 
GaInAS-AlInAS WELL/bARRIER layers CONSTRUCTED TO ACT AS A/4 BrAGG REFLECTORS. 

There is currently much interest in such states in solid state physics [25] . 

9 Conclusion 

We have REVIEWED THE CONCEPTS THAT GAVE RISE TO QUANTUM MECHANICS. THE STATION- 
ary states were introduced first by bohr to explain the stability of atoms and the 
experimental findings on the variation of the electronic current when electrons 
collided with mercury atoms. the generalization of those ideas were discussed 
and the schrodinger equation was introduced. the stationary localized solutions 
of that equation for various potentials with closed and open boundary conditions 
were worked out in detail and the physical meaning was stressed. 

There are other cases of interest, like the solution of the Schroedinger equa- 
tion OF electrons moving in a solid, that are treated in detail in other CHAPTERS. 
Our INTEREST HERE IS TO STRESS THE HISTORIC DEVELOPMENT OF QUANTUM MECHANICS AND 
TO SHOW THE IMPORTANCE OF THE STATIONARY STATE CONCEPT. 
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